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I. INTRODUCTION 


In an earlier work Bohannan [2] developed a mathemati- 
cal model for a system of stored items which are periodically 

- reworked to improve their reliability. That model yields 
an expresSion, for the reliability of an item following 
rework, in the form of a Series. More recently, Bishop [1] 
applied the discrete control theory to systems encountered 
in economics and operations research rather than the more 
typical applications in electrical systems. 

In this thesis the mathematical model for the system of 
stored items which are periodically reworked is developed 
in terms of control theory, and the solution obtained for 
item reliability has a closed form that is more suitable for 
further studies in this area. Important characteristics 
Such as system stability and the steady state value of the 
reliability are directly derived from the control model 
rather than obtained by inspection of simulated values. 

The reliability of the stored items will depend upon 
the reliability at the time of acquisition together with 
the storage environment and length of time the items are 
Stored. When periodic rework of items is done, other varia- 
bles that can affect the reliability or the item are the 
effectiveness of the rework, the rate of rework, the rate 
of acquisition of new items and the rate and policy of 
expenditure of items for use, obsolescence, or even training 


purposes. 








Such a system might include a stock of ordnance which 
1S acquired, stored, and periodically reworked, but not 
expended except for war time use, or alternatively expended 
in training with replacement by new items. 

On the other hand, items such as big missiles are not 
likely to be expended, and thus expenditure and replacement 
cannot be considered as a control variable for the relia- 
bility. Thus, unless due to other factors, the expenditure 
of such items would be avoided in favor of rework of the 
items. 

In this thesis we will study only the case where there 
1s no expenditure and replacement, leaving reliability over 
time to be maintained by rework. We will also consider 
the case where we have several levels of rework and want 
to find the optimal sequence of levels over successive 
reworks in order to optimize a given performance measure. 

A general rework model is developed in Chapter II 
which relates the reliability of an item following rework 
to its reliability following the last rework. An equiva- 
lent model is also developed relating the reliabilities 
of an item immediately before consecutive reworks. These 
models may be solved for any rework cycle given the initial 
reliability at acquisition. 

mai Chapter III the concepts of control theory are 
reviewed and applications to control of the reliability of 


items in inventory are suggested and structured. 





When we have random errors in our meaSure of the item 
reliability or when the rework process introduces random- 
ness in the item reliability, then we have to make our 
judgments based on estimated values. In Chapter IV, a 
Kalman Filter is used with our control model to accomplish 
this estimation. An examination of the final value and the 
stability of our system is followed by an extension to 
more than one dimension for the case where we have systems 
composed of several subsystems in series, parallel ora 
Bombination of both. 

In Chapter V we shall extend our work to the case where 
we can decide among several levels of rework. Using dynamic 
programming, we will show how to derive the optimal rework 
policy for this case in accordance with several performance 
measures, namely, to obtain a desired reliability, to mini- 
mize the costs of several reworks, to minimize the time 
to achieve a desired reliability, or to satisfy a combina- 
tion of all these performance criteria. Conclusions and 


recommendations for further study are offered in Chapter VI. 





Pile eV DEOPMENT Om inh REWORK MODEL 


In this chapter the general rework model is developed 
which will relate either the reliability of an item follow- 
ing a rework or the reliability of an item immediately 
before a rework to the same reliability for the last rework 
and the effectiveness of the rework process. A solution 
is found for these two reliabilities at any instant given 


the initial or acquisition reliability of the item. 


oe THE EFFECTIVENESS OF REWORK 

One form of rework mechanism would raise an item's 
reliability to a certain level which is independent of the 
item's reliability prior to rework, as might be the case 
when components or parts are replaced. Another way, 
developed by Bohannan [2] and adopted in this thesis, is 
for the rework mechanism to achieve an increase in relia- 
bility which 1S proportional to both the item unreliability 
before rework and the effectiveness of the rework mechanism. 
This type of rework mechanism might exist where major assem- 
blies or subassemblies are tested and repaired rather than 
replaced. Denoting the item reliability immediately before 
the rework by R, the reliability following the rework by 


we can define the effectiveness of the rework process 


2. 
lic 
7) 








The interpretation of the rework effectiveness @ is that 
the greater its value, the more effective the rework pro- 
cess or in other words the item unreliability, 1-R, will 
be reduced by an amount proportional to a since from 


Equation (1) we can write 


Under these circumstances, 1t may be possible for such a 
rework mechanism to improve the item reliability to where 
it is “better than new", "as good as new" or “not as good 
as new". Because of deterioration of reliability during 
Storage, we note that in this last case, the reliability 
Might continue to decrease until the items need to be 


replaced, rather than reworked. 


me LHE REWORK MODEL 
From the definition of effectiveness of the rework 
process , given by expression (1), the reliability of 


an item following the rework, Ros is 


29, 
il 


R + a(1-R), 


or 


R = (l-a)R+a. (2) 


10 





Let R(t) represent the reliability of an item of age t 
when there 1S no rework. We assume that all failures are 
random, suggesting (1) that there are no early failures or 
that some form of "burn in" has been used to eliminate 
early age failures, and (i1) that because the items are 

in storage and not in an operating environment, either 
there are no wearout failures or the time to occurrence of 
wearout is much longer than the projected time until the 
next scheduled rework. Thus assuming that the reliability 


R(t) 1S an exponential function, 


q (atbt) 


fees «(. «60 Since the initial reliability Roe On? fe =. 0 


is 


we have 


meee > O and b > 0. The value of the parameter b depends 
On the nature of the item stored and on the storage environ- 
ment . 

If the system contains N items, which are reworked at 


@ constant rate 9, then the period T to "turn over" the 


ab dk 





inventory 1s 


Lo p27 


We begin by looking at new items that are going to be 
reworked for the first time. The duration of the rework 
Process is T/N and the “age” of a new item selected for 
Bework under a First In, First Out (FIFO) policy will thus 
be a =, and its reliability R immediately prior to its 


first rework will be 


If the item is going to be periodically reworked, we 
can rewrite the above equation for the reliability just 
prior to the (k+l)st rework. This will occur at time 


t = (k+1)T-<, and we have 


When k = 0 in this equation, we get 


T 
T ~b(T ~ ¥) 


and from Equation (5) we have then 


I 





Se in other words the initial or acquisition reliability 
may be considered as the "reliability following the one 
rework", or as the initial value for the reliability 
following the rework, R.(k). Evidently there is no rework 
at the time of acquisition and this equivalence is used 
only for purpose of coherence. The reliability just prior 
to the first rework is then only defined for values of k 
equal to or greater than one. 

The process 1S presented schematically in Figure l. 


Assuming that N or the rework rate are relatively large, 


as ~ 
‘ly y~il and 


For simplicity we will denote R(kT) by R(k), and the 
reliability immediately before the (k+l)st rework, R(k+l), 
can then be rewritten as a function of the reliability 


following the kth rework R.(k): 


This expresses the reliability of an item immediately 
before the (k+1)th rework recursively, as an exponential 
Sunction with initial value R, (kT) which is the reliability 
of the item at the beginning of the kth period. 

Rewriting expression (2) for the reliability following 


a rework with this notation, we have 


ibs, 





— eo 


R(k) Be (k) 
FIGURE l. 


Turn over period T 
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R(k+1) R,(k+D 


———— 
Rework Time 


Reliabilities during the kth period 
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Combining the above equations, we have the expression 
mer the reliability after the (k+l)th rework in terms of 


the reliability after the kth rework, 


We may also write the reliability immediately before 
the (k+1)th rework in terms of the reliability just prior 
mo the kth rework, 


nile Tapia. <n ne ee (9) 


R(k+1l) = R(k)e 
When using this equation we have to consider as initial 
Mondition the reliability immediately before the first 
=JoNe 


rework, R(1) = Ro e 


_eeetae SOLUTION TO THE REWORK MODEL 

Equations (8) and (9) are first order difference equa- 
tions, and we can solve them by using the z transform tech- 
Nique, described in Appendix A, with the advantage that we 
Will obtain a closed form solution instead of a solution 
in the form of summation. 

In Appendix A we show that taking the z transform of 


both sides of equation (8), we get 


LS 





-bT Giz 
Ss O S z-tl’ 


where R, (Zz) means the zg transform of R, (k) . 


Rearranged, this gives 





bT ” AZ 
(z - (l-a)e )R, (2) = ZR ~ are i 
or 
Z AZ 
R (Zz) = R ane a. ao Si Se 
: maemo Wi et eae aa 1) 
which when expanded in partial fractions yields 
Z, Z e! Z, 
R_(z) = = ig + — 
7 z-(l-a)e oo oa I Ete.) x= ok [z= 1a) e res 
CL 
Cee) Gro 


The solution of this equation is obtained by taking 
its inverse z transformed, as developed in Appendix A, 


yielding 


@ieniee eat 


a) 
Using the same method for Equation (9), we take the 


meeransform to get 


16 





hg 


-bT 
ipo iieia ie = oR re 
Z-1l 
Bolving for R(z), we get 
Z, ae Pt 
R(z) = = Relea = 
pars Glee eo Eee 


-bT 
mz) = ae (yore or 
eo Slee ica (ee (ele ae 
ee ae eee ae? [(1-a)e 7) 
(Glo ge TI Cas eae eat 


Taking the inverse z transform, we find the solution in the 


time domain to be 


Ql, 


Eguation (13) gives the reliability of an item immediately 


before the rework at time t = kT for a given initial con- 


dition R(1l), the reliability of the item just before the 


first rework. Equation (ll) gives the reliability of the 


item after the rework at time t = kT for a given initial 


condition Ry, the reliability of the item at the time of 


Ly 





SMnisition. Comparing these two equations we notice that 
they are different only by a factor of ees aS might have 
been expected from the Expression (6) that relates these 
two reliabilities. 

As an example, values of R. (x) and R(k) are shown in 
Table (I) for various values of Roe a, b and T. As can be 
seen in Table I, for T = 400 and RS = 0.8 we have, for 
a = 0.8, the case "not as good as new"; for a = 0.569 the 
case "as good as new"; and for a = 0.7, an improvement of 
item reliability constituting the "better than new" case. 
The value of a that leads to the case "as good as new" can 
be found by setting the reliability following the first 
rework R, (1) equal to the initial reliability Ro: Thus 


from equation (11) we have 


RAG) 
or 
a - aR(O)e =  R(0) Geen ee 


sOlving for a yields 
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TABLE I 


Examples of reliability immediately before rework 


R(k) 

ee 761 
0.883 
0.918 


0.928 


0.958 


0.960 


0.961 


0.961 


0.001 for all cases 


ILS, 


R(k) and reliability following rework R,(k). 
immferal reliability = 
Rework Effectiveness = 
Deterioration parameter 
Turnover period = 


R 
fe 
Be 
—1e 
T 
T = 400 
el} 0.7 
\ 
R(k} Re (k) 
0.536 0.861 
OPSSLy 0.873 
0.585 0.876 
0.587 0.876 
= 0.7 
Te —22 50 
R(k) Ro (k) 
0.623 0.887 
O69 0.907 
0.701 0.912 
02710 0.913 
= 100 
Cc — 
R(k) R_ (k) 
0.860 0.986 
0.892 0.989 
0.895 0.990 
0.895 0.990 





DB. THE CONTROL VARIABLE 

As pointed out by Bohannan [2], with these values a 
steady state in reliability is reached by the fourth rework 
and the influence of the initial reliability is lost. Thus 
in order to rework an item to achieve a desired reliability, 
we have to choose values for the rework effectiveness 4, 

Or the turnover period T or both, so that the steady state 
is not below the desired value. The turnover period can 

be shortened by increasing the rework rate 9 or by reducing 
the inventory size N. Another way to have a higher relia- 
bility would be by improving the storage environment, which 
would reduce the reliability deterioration rate b, if this 
is possible. However, the inventory size N and the deteriora- 
tion rate b are not generally considered part of a rework 
Melicy. In our formulation only the rework rate p and the 
rework effectiveness a could be changed within a rework 
BOlicy. 

In this thesis we choose to work with a fixed rework 
rate o and to use the rework effectiveness a to control the 
item reliability. This appears to be the usual practice in 
preventive maintenance involving a fixed schedule and 
several levels of maintenance or rework (corresponding to 


Several values of a). The use of the rework rate op to 
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control the system makes the system non linear, involving 
equations that are time variant. This can be solved by 
another area of the control theory which while not covered 
in this thesis, can be done in a further work. 

In this chapter we derived a rework model for the case 
when we have a constant rework rate. We solved our dis- 
crete recurSive equations for the reliability before the 
rework, R(k), and after the rework, Res) by the z trans- 
form technique. The result iS in a closed form rather than 
in the form of a summation, making this solution more suit- 
able for further analysis as will be done in Chapter IV. 
Before this analysis, we will first derive in the next 
chapter this same closed form solution by applying the 


control theory approach. 
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Iti. DEVELOPMENT OF THE CONTROL MODEL 


In Chapter II we developed the rework model and found 
the system of first order difference equations relating 
the reliability immediately before the rework R(k) to the 


reliability after the rework R, (kK): 


and 


Foie =O (eR KK) ) Rik) s- (2156) 


We then solved this system uSing the z transform technique, 
obtaining an expression for the reliability after any 


Bework as a function of the initial reliability, Roe 


OP? ieee ey 


em, 
and an expression for the reliability immediately before any 
rework as a function of this reliability before the first 


rework: 


Rik) = ({(1-a) 


ee 





In this chapter we will introduce a control theory 
approach to the rework problem, and then derive these same 
expressions in terms of control theory. This will provide 
a basis for later solving the problem of obtaining an 


optimal sequence of rework levels. 


A. THE CONTROL THEORY APPROACH 

In this section we will review, in a simplified form, 
the main aspects of control theorv that are of interest 
to our rework model. 

The first fundamental concept underlying a control 
system is that of a dynamic model, that is a model des- 
Cribing the relationships among the relevant variables and 
parameters of our system. In the model, the variables are 
allowed to change with time in a deterministic way, and 
the system behavior depends not only on their values at any 
instant, but also on their past values and the rate of 
change of these variables with time. The equations derived 
in Chapter II constitutes such a dynamic model of our system. 

The second fundamental concept we will need is that of 
Mee losed-loop control system, i.e., one in which the output 
has a direct effect upon the control action, as shown 
Schematically in Figure 2. The output and the reference 
input are compared in the controller and we say that the 
Output is fed back to the controller. The difference between 
the reference signal and the feedback signal, called the 


Mentor Signal or input Signal to the plant, is then used to 


Ze 
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FIGURE 2. Block Diagram of a Control System 
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change the behavior of our plant or process with the goal 
of bringing the output of the system to a value that gives 
a zero error. 

The "plant" is any physical or abstract object to be 
controlled. Sometimes we feedback a measure of our output: 
for this case, a box representing the measurement device 
would be drawn in the feedback branch of the block diagram 
of Figure 2. In summary, the term "closed loop" implies 
the use of feedback action in order to reduce system error. 
An open loop control system 1S a control system in which 
the output has no effect upon the control action, that is, 
the output 1S neither measured nor fed back for comparison 
With the reference. An advantage of the closed-loop control 
system is that the use of feedback makes the system output 
relatively insensitive to external disturbances and internal 
variations in system parameters, Since a correction action 
is taken anytime the error 1S not zero. 

An useful tool in control theory is a Signal flow graph. 
This is a pictorial representation of a set of simultaneous 
algebraic equations in which each variable is represented 
by a graphical symbol called a node, and the dependencies 
between pairs of variables are represented by directed 
branches drawn between pairs of nodes. These dependencies 
between two variables are called transfer functions, or 
Gains and are defined as the ratio of the incoming variable 


Meetne variable at the end of the branch. Transfer functions 


2D 





are used to label the branches, and indicate that a multi- 
Plication operation is done upon the value entering the 
branch in the arrow direction, delivering a new branch value 
to the node where the branches terminates. The nodes are 
also summing devices which sum all values arriving by the 
way of incoming branches. As an example of this, Figure 3 
shows the signal flow graph for a hypothetical system 


defined by the set of equations 


Xx, = bx. 2 AX ee el 
X, = ax, 

Xx, = ix, + ex, + gu 
ee ae 

Y = hx , 


Given the signal flow graph of a system we often wish 


to write the overall transfer function from an input or 


reference value to the output. 


directly from the signal flow graph by 


This can be obtained 


inspection or 


use of Mason's Theorem. 


This theorem, 


which we will 


by 


use 


later, expresses the transfer function in terms of various 


loop gains and the parallel gains from input to output, and 
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mone 3. Example of a Signal flow graph diagram 


ail 





States that the transfer function from input u to a response 


fy is 
) eae 
Mm = xv = aL 
U A 
where the terms are defined as: 
(1) A is the determinant of the feedback configuration 


Mma is calculated from the equation 
a t oe ’ 
A = 1 .- pL eee ee ees Ete ty 


where the Li (Or Lys etc.) are loop gains all the way 
around a feedback loop in the system. Thus )L. means the 
sum of all loop gains. The next term )' L,L, is the sum 
Meet | products of pairs of different loop gains - e.g., 
LL. and so forth. The prime on the Summation means we use 
meey the products for pairs of gains of non-touching loops. 
In other words, L,L. iS included only if loop 1 does not 
touch loop 3 (two loops touch if they have at least one 
node in common). Likewise, eae e is the sum of all 
products three at a time, where again each of the three 
loops does not touch the other two. 

(2) PB, tomas cache Galion, Ene Ten £OEward path -Erom 
Biput U to output Y (a path which contains no loops). 

f3) A. is the system determinant A after we have 


excluded all loops which touch the PB, Dales 
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The Mason's Theorem is better understood by an example. 
Suppose our system is described by the signal flow graph 
of Figure 3 and we want to determine the overall transfer 


function Y/U. The loop gains are 


L = ab, L = ce , if = alicd 


Hence 


)L. ab + ce + aicd 


The Loop L. touches both Ls and Lo Due L, and L. are 


non-touching. Hence 


-, oe 
eb = aided, 


and there are no three non-touching loops. We can then 


Mmeaite the determinant A 


[> 
| 
4 
i 
eat 
a, 

+ 
Bat | 
is 
aw 
i 
rd 


Or 


A = 1- (ab + ce + aicd) + aicd. 


mie Girect paths from U to Y¥ are 


5) 





and 


Since Py touches 


Path P~ does not 


2 
2 
SuDStituting 
At 
ex 
At 


ee ten ae aia) 


= oe it 


all loops, 


EOuch 1 hence 


le’ 


into Mason's formula, we have 


Pi4, x BoA. 


A f 


Ffaich + ch(l- ab) 
1 - (ab + ce + aicd) + aicd 


Summarizing, we can say that Mason's theorem gives a simple, 


Fast procedure for writing a desired transfer function 


directly from the signal flow diagram. The response of 


Bur system can then be immediately derived from this trans- 


mee function as proportional to the input, 
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and this is called the solution to our control model, 
i.e., an expression for the response of our system as a 
function of any input to the system and a known (derived by 


Mason's Rule) overall transfer function. 


me CHE RELIABILITY CONTROL MODEL 

Our concept of a system is not limited to physical 
systems. The concept can be applied to abstract, dynamic 
phenomena such as those encountered in economics and opera- 
tions research. Also, feedback control systems are not 
Mmimited to the field of engineering but can be of particu- 
lar interest to the manager, public official, operations 
researcher, biologist, and design engineer. In this section 
we will structure the reliability model in control theory 
terms. The plant of the control theory approach will be 
the inventory, where the (decreasing) reliability of items 
is our variable of interest, or state variable. The con- 
troller will be substituted by a decision-making function 
that will decide upon a level of maintenance which is optimal 
in the sense of minimizing a certain performance measure. 
Memewill defer optimization until Chapter VI, and for the 
present will assume a single level of rework, as defined 
my Section D of Chapter II, corresponding to a certain 


effectiveness of the rework a. 


Sy 1h 





Our dynamic model is described by expressions (15) and 


(16), 


and 


and the signal flow graph can be drawn in several ways, one 
of which is shown in Figure 4. At the second node from 

the left the present item reliability, whose measured value 
is fed back to this point, is compared with the reference 


value, R that is made here equal to 1.0 due to our 


D! 
definition of rework effectiveness (this will be better 
Clarified later). The difference, when not zero, produces 
an automatic decision to rework the item, and after the 


rework, this difference becomes multiplied by a since from 


expresSion (16) we have 


This difference will be denoted by U(k), and called 
the input to our plant. Physically the item is returned 
to the inventory at this point and the reliability of the 


item is given by the sum 


SZ 
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FIGURE 4. Reliability Flow Graph for the Rework Model 
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This sum is performed at the fourth node. While in the 
inventory the item reliability will deteriorate by a 


oT given by expression (15) and at the end of a 


factor e. 
meeod Of length T we have the item reliability R(k+1) 
immediately after the (k+l)th rework represented at the 
fifth node by R(k+l). The relation between R(k) and R(k+1l) 
is a function called a unit delay of time, represented by 
a box in the flow graph, with the meaning that the output 
of our system R(k+1) 1s one period T ahead of the value 
Ged back, R(k), used to compute it. 

The variable Y(k) is used to emphasize the fact that 
R(k+1l) is the output of our control system. Branches with 
a gain equal to 1.0 are used when there are variables that are 
just renamed or going to be operated on in another mode. 

Taking the z transform on both sides of Equations (15) 
and (16) by applying the properties given in Appendix A, 


we have 
R(z) = z R. (2) ieee ape in (20) 


and 





and the signal flow graph can be drawn as in Figure 5. 
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BOGURE 5. 


y Ulz) Ll Rez) eT zR(z) zl JRL) 1 


Reliability Signal Flow Graph to the Control 
Model in the z Domain 


5 





The response of the system, Y(z), can be now obtained 
by applying Mason's theorem as shown earlier. Since we 
can consider that we have here two input values, namely 
R(1) and Ry: this yields the overall transfer function 


with respect to the initial condition R(l1) as 


Nez) ak Z 


Rees ees & eel peaen a 





and with respect to the reference value Ry: 


= DPT o.-Z 


ae ( Te 


ome ee 


Rp ae ee 








(i) Poe eens) 


The response of the system due to the reference value 


that here is made equal to 1.0 is then 


=65 7 
R52e Z 


ee en ome 


) 
and the response of the system due to the initial condition 


is given by 


wie) = REE) 2 _ 


(z-(l-a)e ~~) 


Using the superposition principle for differential 


equations, the total response is 


" ae oe Rielz 
7 a > == 
Ae eeega ESO) ) 2a eee ) 
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This is the same as the z transform of Equation (13) 
in Chapter II, and the inverse z transform of the total 


response then gives the same Equation (13) for the relia- 


bility of an item immediately before a rework. This is 
k -bkT ae Pt k _-bkT 
Y(k) = R(k) = (l-a)e Ras) pr [1l-(l-a)e 
i= )e Bek 
k = | pee ae ee A 
Where the argument (k+l) was changed to k. In this case 


the z transform had to be applied because the unit delay 
present in the system, does not provide a linear function 
necessary for the application of Mason's theorem. In the 
z domain, however, the unit delay becomes a very simple 
function and the solution may be found easily. 

If the reference input Ry cS not cequale to 1.0 bites 
the desired value of the item reliability that we want 
R. (k) to match after some reworks then the definition of 


mework effectiveness would become 


Here, if the reliability immediately before a rework was 
greater than the value of the desired reliability, the 
Value of a would be negative. To avoid this problem we 
Might use a decision rule that could be Ee nessed asa 


eunction M(R, R,) defined as 


a7 





with the meaning that no rework has to be done when the 
item reliability R(k) 1S greater than the desired value 


R.. Our model would then become 


This approach, however, would introduce a non-linearity 
into our system that could not be eliminated by the z 
transform procedure, and the solution of our problem would 
probably have to be obtained by computer. Hopefully we 
can handle this problem and find an optimal sequence of 
reworks that makes the value of Re (x) to match the desired 
value Ry by uSing the optimization process that will be 
developed in Chapter V. 

Summarizing, we have structured the item reliability 
Model in terms of control theory for periodic reworks with 
effectiveness of rework a as defined by equation (1). We 
are now able to apply the tools of control theory to find: 
(1) the limits of stability of our system as a function 
Of the rework period T, (11) the stable final value of the 
reliability as the number of reworks increase, and the use 


Of a Kalman filter. This will be shown in the next chapter. 
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ne oe Rebeca TONS (Or THE CONTROL THEORY MODEL 


In this chapter we will show how the results of our 
control model can be applied to study the behavior of the 
item reliability. First, we will study the stability of 
Our system, and then the final value reliability when t 
increases. After introducing randomness to the system, 
we will show the use of a Kalman Filter to estimate the 


value of the reliability that is fed back. 


me oOYSTEM STABILITY 

A linear control system is stable if the output even- 
tually comes back to an equilibrium state when the system 
is subjected to a disturbance. In our reliability system 
we are interested in knowing if there 1s any range of values 
of the rework effectiveness factor a for which the system 
becomes unstable, i.e., the rellability might continue to 
decrease until the items need to be replaced rather than 
reworked. 

From control theory we know that the characteristic 
equation of a control system is equal to the denominator of 
the overall transfer function plus 1.0 and that for the 
system to be stable the roots of the characteristic 
equation, in the z domain, have to be inside the unity circle 
[10]. For our reliability system we have from Equation (21) 


the characteristic equation in the z domain as 
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AS1D) Tees Glee aes eel ee Os 


Rearranging our characteristic equation, we have 


a= Cloner ne 1s isan ome an eu 


With complex roots at 


-bT 
y= tie + 4 5 Be) eee locate 


~bT, 


The absolute value of these complex roots 1s given by 


[z| = S[(1-a) “2 *P7-1} 
Here, 

jz[ < 1 
or 

jzZj° <1 
when 

=[ (1a) “e a I ae 
or 
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that when solved shows us that for a stable system, the 


value of a must be within the range 


For the example of Chapter II with b = 0.001 and 


T = 100, we get 


Oe 36) 67a F<) 22515 


that agrees well with the value of a in Table I for which 
the item is "as good as new" after each rework. 

Outside these limits or when a is less than 0.56, (since 
a®@ by definition is less than unity), the system is unstable. 
Here, the item will have its reliability deteriorating with 
time until the item needs to be replaced rather than 


reworked, li.e., for a < 0.56 eventually R- 0 


Bm THE FINAL VALUE IN STABLE SYSTEMS 
When the system is stable its response tends to a 


constant value as the time increase whenever there 1S no 
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disturbance input. To find this steady state value we 
will apply the final value theorem, derived in Appendix A, 
that gives us the limit of the system response when time 


becomes very large. The theorem is stated as 


iam Y tk) = lam [(z2=<1) Y(z) } 
k>+ Z>l 


For our overall transfer function, this results in 


GZ ZR(0) (z-1) 


iim Ro UK) Norte 
ie > Z>l 2-(l-a)e oe z-~(l-a)e oe 
O 
= - (26) 
ei anew ss 


Substituting selected values from Table l, we get 


a = 0.7 T = 100 jam Bek) ee= WOL96L 
tC7o = 

m= 0.8 T= 00 Lim R. (k) = Leeroy 7, 
t-+c 

a = 0.7 T = 400 lim Ro (k) = 0.876 
C>c 


and these limits agree well with the values taken from 
the table for the fourth rework. 
As we can clearly conclude from expression (26) the 


mecady state value of the reliability after the rework R.() 
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depends only on the rework effectiveness qa and on the turn- 
over period T for a given degradation parameter b, and is 
thus independent of the initial value of the reliability 


as waS intuitively observed in Chapter II. 


C. ERROR COMPENSATION USING A KALMAN FILTER 

In this section we will show how a Kalman Filter can 
be applied to our model in order to compensate for errors 
in the measurement of the item reliability. Again, measure- 
Ment here does not imply a physical action. We suppose the 
nature of the stored item 1s such that it receives frequent 
diagnostic checks which yield an estimate of its reliability, 
and such "measurements" are subject to random error. 

In the model developed in Chapter III we have shown how 
the item reliability is fed back so that a decision can be 
Made about doing a rework so that the item's reliability 
can approach a desired value. In practice we have to 
measure the item reliability or estimate it by some way in 
order to feed back this observed value. This measurement 
process, however, contributes to the variance within the 
System. Other sources of variance in the process include 
Variability in the item environment and in subsystems 
reliabilities due to repairing or substitution of parts. 

We call these random errors as noisy or random input to 
Our plant, since we can model these errors as random 
Variables that are input somewhere in the plant as will 


See later. 
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If we do not want our control system following these 
random inputs excessively it would be a good policy to 
smooth the reliability values over time and then feed back 
a predicted value upon which we made our decisions. 

A useful, smoothing device is a Kalman Filter, and a 
block diagram of a general system with a Kalman filter is 
shown in Figure 6. Other filters could be used here but 
with a discrete system we chose the Kalman filter: it is 
the optimum recursive filter in the sense that it minimizes 
the variance of the estimator error. The random input at 
the plant in Figure 6 is designated by the random variable 
w({k) and represents the error introduced by the environment 
and rework process. The random error in the measurement 
is designated by the random variable v(k). 

Making the assumption that these noises are additive 
expected values, we can restate our system equations (15) 


end (16) as 


and 


R_(k) = (l-a) R(k) + aR,» 


and substituting the second equation into the first, we 


have 
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FIGURE 6. Block diagram of a general control 
system with Kalman Filter 
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R(k+1) = eT [ (1-0) R(k) taRy] + w(k) 


Rearranging and grouping terms we have 


where 
Des = afR, SR e)) 


as before and Ry = alive 
If we introduce the Kalman Filter in our system we 
are feeding back an estiamted value R(k) instead of the 


actual value R(k) and the input to the plant becomes 


where H can be any function describing the measurement, and 
v(k) is the random variable in the measurement process. 
We need the following additional assumptions to make 
Ehe application of a Kalman Filter valid: 
(a) The measurement error v(k) and the random process 
input w(k) has zero mean and are uncorrelated between 


States, with variances S and Q respectively, or 
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Pe aeonyadel — SOs, k,)] = 0,1,2, «.. , 


E{w(k)w(j)] = Q 6 fo = (Ogl 27 eee iy 


kj ’ 
where 9165 icvtenowkronecker delta fLunction .- 

(b) The initial state is a random variable which has 
known mean E[R(0)] = Ry and variance M. 


(c) The estimator is characterized by the linear 


relationship 
R(k{k) = R(k{k-1) + G(k) [Y(k) - HR(k|k-1)] , k = 0,1,2,.. 
where: 


R(k | Kk) is the optimal estimate of R(k) given 
observations at times up to and including k, 
R(k|k-1) is the optimal prediction of R(k) given 
observations at times up to and including k-l, and 
G(k) 1s the gain of the Kalman Filter at each step. 
(d) The random input and initial state are uncorrelated. 
(e) The random errors in the plant process or random 
input w(k) and the random error in the measurement 


process v(k) are uncorrelated or 


Bey (a) i =O j,k = O,1,2, 


47 





The functioning of the Kalman filter is by predicting 
the value of the item reliability R(k{k-1) using the esti- 
mated value at the end of the (k-1)th period, multiplying 
this predicted value by the function H, so that this new 
value can be directly compared with the measured value. 
Bne filter gets the correction term [Y (k) -HR(k|k-1) ] 
that can be conveniently weighted by the gain G(k) to 
correct the predicted value, and obtains a new estimated 
Malue for the kth period. 

Notice that if we did not use a predicted value, 
what is equivalent to say that R(k|k-1) = R(k-1/k-1) = R(k-1), 
we would have exponential smoothing and the gain G(k) 
would then be the smoothing constant. With the predicted 
value any trend 1s recursively incorporated in the esti- 
mated value. 

The best prediction R(k|k-1) that can be made is 
clearly by using our model Giese (2); 


be 


R(k[k-1) = e@ 7 R(k-1|k-1) + e >? 


Ds h)e, (30) 


and from assumption (c) above, the equation for the esti- 


mated value R(k|k) is 


R(k) = R(k{[k) = R(k{[k-1) + G(k) [(Y¥(k)-HR(k[k-1)]. 
(31) 
With these two equations and the assumptions above it 


is shown [9] that the sequence of gains that minimizes the 
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variance of the estimator error P(k|k) is given by the 


recursive equations 


G(k) = P(k|k-1)H[H*P(k|k-1)+S] > ; (32) 
P(k|k) = [1 - HG(k)] P(k|k-1l) , (33) 
and 
Dalia Sem (eu yr (cel k=) eo: (34) 
These are initialized by the value 
P(O|-1) = M. (35) 


From these equations we can see that the Kalman Filter 
gain G(k) does not depend on the measurement values Y(k), 
and can thus be computed in advance and stored for later use 
when processing the measurements as they become available. 
We can see also that the gain at time k, G(k), is inversely 
realted to the variance of measurement error S — the more 
uncertainty in the measurements, as reflected by a larger 
S, the smaller G(k) will be, and the less the correction 
term fy (k)-HR(k|k-1)] in Equation (31) will be weighted in 
determining the next estimate. The random input error Q 


also affects the gain G(k), and as the uncertainty in the 
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prediction of P(k|k-1) increases and the more uncertainty 
we have in our model, the larger will be G(k). For large 
values of M or no confidence in the initial guess R(0|-1) 
the filter makes R(0|0) equals to the first measurement 
Y(0) since G(0) is equal to 1.0. When the initial guess is 
hetter the filter makes a weighted combination of the guess 
R(0|-1) and the first measurement Y(0), as one can conclude 
from Equation (31). 

A flow graph in the time domain of our system is shown 
in Figure 7. The lower part of the flow graph is the 
Kalman Filter, which receives as input the measured value 
Y(k) and the input to the plant U(k), and has as output the 
estimated value R(k) . This 1S multiplied by the measurement 
function H and the result is compared with the observed 
value Y(k). The difference is then multiplied by the gain 
G(k) and summed to the predicted value to produce the 
estimated value, as given by Equation (31). 

When we have more than one state variable, an advantage 
of the Kalman filter not shown here is that we do not need 
to measure all state variables but only a smaller number 
to get all the estimated values we need to feed back to the 
plant. 

As an example, suppose we have values from Table I in 
Biat © = 0.7, b = 0.001, T = 100 and Ro = 0.8; suppose also 
that the variance in the estimated value of Ro is M= 0.1, 
and that the variances in the random input and in the measure- 
ment are known to be Q = 0.1 and S = 0.1. Let our measure- 


ment equation be given by 
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PIGURE 7. Flow graph in the time domain for the 
rework system with a Kalman Filter 
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Y (k) 


which in terms of Equation (29), 


R(k) 


+ v(k) 


values we can compute the Kalman filter gain for four 


periods using Equations 


G (0) 


C1) 


S23} 


G (4) 


II 


2) ete Omen!) 


Suppose the measured values 


GL) 


U2) 


ale 


, and 


and we guess a value for the initial reliability 


mor |/0) = 0.724. We obtain, by using Equations (30) and 


the values 
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implies H = 1. With these 


(aie, 





Ree  O.7 35/7 ; 


2 ee Ole 2-5), 
Po) pe de oo on ane 
R(4) = 0.856 


Notice that the reliability before the rework, R(k), 
is not defined for k = 0. For this reason we use as initial 
guess R(1]0) the value of the reliability before the first 
rework computed from the initial value Ry = 0.8 and use the 
gain G(0) to estimate R(1), the gain G(l) to estimate R(2) 
and so on. The gain would not be delayed by one period if 
we had decided to call the reliability before the first 


rework by R(0). 


D. EXTENSIONS TO MORE THAN ONE DIMENSION 

In a missile we can consider for example its booster, 
Cruise motor, guidance, and warhead as four subsystems that 
are independently influenced by the environment and have 
different reliability deterioration factor b. To extend 
Our model for this case we need to decompose our system 
in several distinct subsystems. Here we may put the 
equations for the subsystems in matrix form. 

We can define the vector of reliabilities of the 


Subsystems that will be called our state variables, as 


De) 





meumeeel Vector of reference reliabilities as 


| The input vedtor can now be defined as before in terms 


Of the effectiveness of the rework hs for each subsystem: 


04 (k) oe a 
| | 

U. (k) _ R 
pxk) = ; = 10, G5 3 %4! a : 

| U2 Ck) Sao 
0 6k) Sy 

. 
(ik) =o LR, en) 
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We can then write our state equation for the system with 


random input as 
R(k+l) = ® R(k) + 9 WK) TS, (36) 


where the matrix @¢ 1s 


=b. 0 
e 1 0 0 0 
-b.T 
0 e 7 9 0 
aa oe 
0 0 2 0 
Jo) ab 
0 0 0 Bie 
a 
The measurement equation is 
ey > HORE) eye), (37) 
and the vector of estimated values is given by 
Comes RK Ky RR Red) tC (kv (Rye = RCE) 
(38) 
and the vector of predicted values by 
R(k|k-1) = e@ 2° R(k-1[k-1) + @ 2” U(k-1), 
(39) 


Where G(k) is the vector of the Kalman filter gains. 


Bie 





If we are interested in the overall system reliability 
and if as in the missile case the subsystems are inde- 


pendent and in series, we have the output equation 


which can be added to the control signal flow graph as 
shown in Figure 8. For subsystems arranged in parallel 

or series-parallel, the system can be handled in a similar 
manner, and needed changes in the output equation will not 
influence the system itself. 

In this chapter we applied tools from control theory 
to our reliability model and showed how to find the limits 
of the rework effectiveness a for which the system is 
Stable and the final value of the item reliability when k 
becomes large. We also showed how a Kalman filter may be 
applied to smooth the estimated values of item reliability 
and feed back this estimated valuesupon which the decision 
about a 1S made. 

In the next chapter we will show how the decision 
about a can be optimized with respect to a given performance 


moaSure. 
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FIGURE 8. Signal flow graph for multidimensional rework 
control system 
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V. OPTIMIZATION 


In the previous chapter we have studied our control 
model for the case we have linear equations and a single 
level of rework effectiveness a. In this chapter we will 
extend this study for non-linear equations and several 
levels of rework effectiveness, and show how to find the 
sequence of rework that optimizes a given performance 
criteria. First we will discuss various performance 
Criteria that could be used for our control model, and then 
we will show how to find the optimal sequence of rework 


levels by dynamic programming. 


A. THE PERFORMANCE MEASURE 
The optimal control problem in our thesis is to find 


a control u(k) as defined by the expression 


u(k) = a(k) (Ry - R(K)) , (40) 


which causes the system to follow a trajectory that opti- 
mizes a performance measure J. The rework effectiveness, 
a(k), is now allowed to change from step to step and in 
fact the problem now becomes one of finding the sequence 
of rework effectiveness a(k) that yields the optimal 
Sequence of control u(k). NOtice that since a(k) is no 


longer constant, Expression (40) is non-linear. 
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First let usS consider a trajectory (that is, the 
sequence of values the item reliability follows over time) 
that attains the desired reliability with a smallest num- 
ber of reworks. This type of problem is called a “Minimum- 
Time" problem andthe performance measure to be minimized 


may be generally stated as 


for the case of our discrete system, this becomes 


Wewhere (N-1) is the first period the desired reliability, 


Ry 


is attained. 

If instead we want to minimize the deviation of the 

first state of our system from its desired value we have 
the type of problem called a "Terminal Control" problem. 


Here, possible performance measures are 


N 
a tee ee OS) 7 


ae 


where the summation is done in all dimensions or reliabili- 


ties of the sub-systems. If positive and negative deviations 





3), 





are equally undesirable these deviations should be 
squared. 
To transfer a system from an arbitrary initial state 
Ry to some specified desired value with minimum expenditure 


of control effort or with minimum cost, we need to minimize, 


N-1 N-1 
ar eda) r =n) souds)r 7, 
k=0 k=0 


where we assume that the cost of a rework C(K) 1S propor- 
tional to the control u(k) or equivalently that the cost 


ms proportional to the reliability improvement, R_(k) - R(k), 


Ss 
| after each rework. 

Since these criteriaare completely distinct in their 
: concepts, we Shall follow the usual approach of using a 


' combination of them, that for the one dimensional problems 


takes the form, 


N-1 
SS sa) eR se (1) ae te eee dat 
k=0 
Mor 
N-1 
J = HN[R, pe) Y [Q + Gu(k)] , (41) 
k=0 
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where HN, Q and G are relative weights so that by adjusting 
their values we can weight the relative importance of each 
criteria with respect to the other. These weighting fac- 
tors can be functions of time if the relative importance 
varies with time. Notice that now the performance measure 
is a combination of criteria and has no physical meaning. 
It should be noticed also that the rework effectiveness 
a(k) now has the definition given at the end of Chapter III: 


Roe) = Ri) 
COS) ee ee et 


Re = Rie) 
and that even when Ry 1S not equal to 1.0 we do not need 


the deciSion rule 


because the constraints in a(k) do not permit a to become 


negative. 


B. OPTIMIZATION USING DYNAMIC PROGRAMMING 

We now wish to show how the performance of our rework 
control system may be optimized in terms of the combined 
performance measure (41). An optimal solution would be 
expressed as a sequence of reworks that minimizes this 
Criterion. We shall approach the problem with the method 


Of dynamic programming, as developed by Bellman. 
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The basic notion here is given by Bellman's principle 
Bf optimality: 

"An optimal policy has the property that whatever the 
initial state and initial decision are, the remaining 


decisions must constitute an optimal policy with regard 


to the state resulting from the first decision" [5]. Thus 
af Seay ny (R(KFL) ) is the minimum cost to reach the final 
4 


State at k = N, starting from the state R(k+l) at time 

t = (k+1)T, and Sia teat is the cost to move from kth to the 
(k+1l)th state, then Oey (ROK) ,u(k)), the minimum cost to 

go from the kth to the (k+1)th state when we use a particu- 


lar control u(k) 1s given by 


ce (R(K) ,u(k)) = J (R(k) ,u(k)) + 


ical CEC )0)a2 


* 
Jel, N 
(42) 
The optimal decision at the instant k, u*(k), is the 


decision that minimizes Cy (RIK) pu (Kk) ) over the set of 


possible controls u(k), 


Jy (RIK) ) = min [C¥ (RUK) -u(k) ) J ‘ (43) 
u (k) 
Expressions (42) and (43) form the functional equations 
Of our dynamic programming approach. The optimal sequence 
of decisions are built up from the final state N backwards 
toward the earlier states. This is necessary in order 
that J* be known prior to the calculation of C*. The values 


of R(k) are given by our model, and the constraints of the 
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problem are 


© 
© 
A 


R(k) 


[A 
- 
S 


and 


The first step in the computational procedure, then, 
is to find the optimal policy for the last stage of opera- 
tion. This is essentially a matter of trying all of the 

allowable control values at each of the allowable state 
values. To limit the required number of calculations, and 
| Make the computational procedure feasible, the allowable 
State and control values are discretized. The degree of 
approximation depends on the separation of these discrete 
values and on the method of interpolation used and can, 
of course, be adjusted. 

For each discrete value of R(N-1) we try all discrete 
' values of u(N-l) and calculate the resulting state R(N). 
The optimal control for this rework is the one which yields 
the minimum cost. The procedure is repeated for all the 
Other discrete values of R(N-l). This gives a table of 
Optimal policy for each value of R(N-1) at the last stage. 
Since the cost J 1s dependent on the value of the 


N-1,N 
State R(N-1l) and on the value of the input applied, u(N-l), 
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the minimum cost J* and the optimal control u*(N-1) are 
dependent on the value of the state R(N-1). For the last 
rework we use only the first term (terminal control) of 

our performance measure (41), but for the other, successive 
stages we must compute each term in that expression. When 

a state does not coOlncide with one of the discrete values 

we have to use interpolation to find the corresponding value 
of the performance measure. Since a direct search is used 
to solve the functional equations, the solution obtained 

is regarded as the global minimum. 

Petlow chart describing the computational procedure is 
maown in Figure 9. A Fortran program for this flow chart 
is described in Appendix B. 

As an example of the use of this dynamic programming 
procedure, we sought the optimal sequence of reworks for 


the case where b = 0.001, T 


100 and. constraints; 


One <> Roe <0 


and 


Gee on atk) i eae 


We chose to work with twenty discrete values for R(k), 


thirty discrete values for a(k), and four reworks. We 


input in the program the equations of our model: 
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; READ input data: 

I staces = = N; number of state 
lwalues = NX; number of 
‘control values = NU; other | 
recuired information j 


| 
| 


number of 


X= 9 
. at 
Caiculate and store Jia: | 
(minimum cost of last stage) 
£5r ail discrete valimas of ' 


Rts) | 







x+1 
Sat Rk) 5 equal *o the lesser 
ciscretizec value by ms 
'y — 


H:? 


Sasicine Mawer: Son U(s) 
equal to the lesser 


| Set COSMIN to a large 


Beet ee ue by 
1 





CALCULATI the value of R(N-xX+1). 
free =<S value to select the | 
| avorcoriate stored value of | 
| a) iN (enters er Securzed | Caance U(K 
Mee 190 COlincides with a discrete , 

i 

] 





| 
‘ 
hem See = 
; fee) | CO Soa Se rien sost discre+=2 
j Over she final (eX) stages for he = 748) 
i UiN=X) | 









ss 


COSTOP (N-K,3) and 


. ' hae 
\ } = “ a6 
te a } ; EN sty) UCE - (N= t‘~ ) 


Gocre Coola, 11h 
< 






(“% 


Pe 
- — 
a 


“jess “nan 
COsiiN 2? 









Yes 


7 
: 


STORE the valve C* in 
COSMIN; STORE che value 
ign 27) ‘MIN 





ere ee eed es 


COSTo2 (N+ ? a sor 
lail discret points! 


————— pand all staces. : 
1 ry i 


FIGURE 9. Flow chart of the computational procedure we)ig 
the dynamic programming solution for optimal 
control of the rework control system 
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R.(k) = R(k) + a(k) (1-R(k)) 


R(k+1) = Re (k) e , 


and considered as input variable the rework effectiveness 
a(k) Since in this program the equations need not be linear 
and this is the variable of immediate interest for the user. 
For the performance measure J given by Expression (41) 
we considered four distinct types of problems obtained by 
changing the weights. These are: 

(1) Minimum Cost, so that the weights of the combined 


Criterion J (41) are 


Cau. 0 ; Q = 0.0 and HN = 0.0 , 
(2) Terminal State, or minimizing deviations from the 
final value. Here the weights in J are 
G = 0.0 , QO = 0.0 and UN nq can Oey 
(3) Terminal State and Minimum Cost, where we assume 
a desired reliability Ry = 0.99, and a relatively important 


terminal state represented by the weights 
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and 
(4) The linear combination of the three criteria 


represented by the weights 


G= 90.05 , Q = 0.05 and HN = 500.0. 

The dynamic programming solutions for the four cases 
are presented for each rework period in the Computer Output 
Section of this thesis. From these solutions we can now 
find, for each of the four cases, the optimal sequence of 
reworks for a given initial reliability Ro (entering with 
this value at the first rework table and, interpolating 
among the vaiues listed, we can follow the process until 
the final state is reached). 

The results are in Tables II to V for two values of 
initial reliability of an item. As we can see from these 
tables, the sequence of reworks that minimizes our perfor- 
mance measure depends on the value of the initial reliability 
Ros for given values of the desired reliability RB and of 
weightings in the performance measure. 

Summarizing, we have applied a dynamic programming 
procedure to our control model of the reliability of an 
item and were able, for the four selected cases, to find 
an optimal sequence of reworks that minimizes a selected 
performance measure. For other selections of performance 
Measure J, relative weightings HN, Q and G, constraints, 


initial reliability Ro and model parameters b and T, we 
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TABLE Il 


RELIABILITY FOLLOWING REWORK FOR INITIAL RELIABILITY 
Roe AND OPTIMAL SEQUENCE OF REWORK LEVELS a (k) 


Case (1) - Minimum Cost 
Net 
PERFORMANCE MEASURE: J = [R,-R,(k)]H + ) Q+GU (k) 
k=0 


with H= 0.0, Q2+=0.0, G=1.0 and U(k) = a(k) 


R, = 0.8 
k R. (x) ck) Ba Garde) 
0 0.800 0.60 0.890 
1 0.890 0.60 0.924 
2 0.924 0.60 0.933 
3 0.933 0.60 0.936 
4 0.926 
R, = 0.9 
k R_ (Kk) a (kK) Rq (k+1) 
0 0.900 0.60 0.926 
it 0.926 0.60 0.933 
2 0.933 0.60 0.936 
3 0.936 0.60 0.937 
4 0.937 
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TABLE SIit 


RELIABILITY FOLLOWING REWORK FOR INITIAL 
RELIABILITY R_, AND OPTIMAL SEQUENCE OF 
REWORK LEVELS a (k) 


Case (2) - Terminal State 
N-1 
PERFORMANCE MEASURE: J = [R,-R,(k)]H + ) Q+GU(k) 
k=0 
with R, = 0.99, H=1.0, Q=0.0, G= 0.0 


and U(k) = a(k) 


= 0.8 

R, (k) a (k) Rg (k+l) 
0.800 0270 OR 2 

Oo 2 Oe: O2983 
0.935 0290 O39 0 
0-950 0290 02990 
O2990 

Ro = 0.9 

Re (x) a (kK) R, (k+l) 
O00 0.90 0.982 
0.982 0.90 O29 39 
0.989 ONS 0 Or 0 
0.990 0220 0.990 

0 S20 
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TABLE IV 


RELIABILITY FOLLOWING REWORK FOR INITIAL RELIABILITY 
Ro AND OPTIMAL SEQUENCE OF REWORK LEVELS a(k) 


Case (3) - Terminal State with Minimum Cost 
N-1 

PERFORMANCE MEASURE: J = [R,-R,(k)]H + ) Q+GU (k) 
k=0 


with Ry Veo, “2t— 50050 >, “OC = 0-075, 7G = 0-9 


and U(k) = a(k) 


Ry = 0.8 
k R, (k) oS) Re (k+1) 
0 0.8 0.60 0.890 
Ll 0.890 0.60 O2924 
Di 0.924 0.60 0.933 
3 0.933 0.867 0.979 
4 0.979 

R = 0.9 

O 
k Re (kK) a (k) ccm) 
0 0.900 0.60 0.926 
Ll 0.926 0.60 0.932 
2 0.933 0.60 0.936 
3 0.936 0.864 0.978 
4 0.978 
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TABLE V 


RELIABILITY FOLLOWING REWORK FOR INITIAL RELIABILITY 


R,, AND OPTIMAL SEQUENCE OF REWORK LEVELS a (kK) 
Case (4) - Terminal State in Minimum Time with 
Minimal Cost 
N-1 
PERFORMANCE MEASURE: J = [R,-R,(k)]H + ) Q+GU(k) , 
_ k=0 
with Re 0299" SHe= 500.0 7 © = 0205 .. G = 0205 
and Ut) = Gtk) 

R =i 50) a8 

O 
k R, (k) a (k) Reece) 
0 0.800 0.60 0.890 
1 0.890 0.60 0.924 
2 0.924 0.868 0.977 
3 O27 7 0.90 0.989 
4 0.989 

R =. Oo 

O 
k R, (Kk) a. (k) Rq (k+1) 
0 0.900 0.60 0.926 
1 0.926 0.60 0.933 
2 0.933 0.835 0.973 
3 0.973 0.90 0.988 
4 0.988 


(ek 





would obtain different sequences of rework effectiveness 
a(k) that optimizes the selected performance measure. A 
Suitable selection of the performance measure and relative 
weightings being a responsibility of the manager according 
to his purposes, experience and data available. Thus the 
cost of each rework, the minimum value allowable to the 
item reliability, minimum and maximum levels of rework that 
are physically feasible to implement, and other factors 
like the urgency to attain a desired reliability in minimum 
time, Or a constraint in the budget, are taken into account 
when selecting the weightings in the combination of all 
the criteria. 

In the next chapter we will present the general éonclu- 


Sions of this thesis and suggest the areas where one could 


do further study involving the rework control model. 
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VI. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 


A simple, closed-loop inventory storage system from 
which the stored items are removed, sent through a rework 
mechanism to improve their reliability, and returned to 
Storage has been investigated in this thesis. Such a sys- 
tem might be a stock of large ordnance which is acquired, 
stored and reworked at regular intervals. 

A general rework model was developed by using the tools 
of control theory and a solution in closed form was found 
that permits further analysis of the system in a much eaSier 
way than if we had a solution in the form of series. This 
model also permitted us to find in a direct way the values 
of effectiveness of the rework that makes the system unstable 
and the final value of the item reliability over several 
reworks when the system is Stable. 

A Kalman filter was used in the control model to obtain 
an eStimation of the item reliability when we have random 
error in our measure of the item reliability or when the 
rework process introduces randomness in the item relia- 
bility. This model was extended to more than one dimension 
for the case where we have systems composed of several 
Subsystems in series, parallel, or both. Finally a study 
was made for the case where we can have several levels of 
rework and want to find the optimal sequence of such reworks 


that minimizes a performance criteria that leads to the 
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case of minimum time, minimum cost, terminal State ora 
linear combination of these criteria. 

This thesis shows that control theory may be applied 
for administrative problems where a mathematical model des- 
cribing the system can be derived. Control theory gives 
a broader understanding of the problem and due to its flexi- 
bility would permit us to include in the model many other 
variables not studied in this thesis. 

Among the several interesting areas which might be 
pursued in further study, one is to determine methods for 
measuring or estimating the item reliability by frequent 
diagnostic checks so that this measurement process can be 
incorporated in the model. Another area iS to broaden the 
field of the rework model so that other state and input 
Variables could be incorporated. These include the proba- 
bility that an item chosen at random might have reliability 
exceeding some predetermined reliability requirement, and 
Other input variables (Such as items expended for training 
purposes, use, or obsolescence, as well aS new items acquired 
periodicaly). The rate of rework, the rate of expenditure, 
Or the rate of acquisition could then be sought as control 
Variables besides the rework effectiveness. Non-linear 
functions caused by the decision rule of doing a rework 
Only when the item reliability is smaller than the desired 
Value, can also be solved in terms of non-linear control 


theory. 
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A reliability control model has been developed in this 
thesis for a system of stored items requiring rework. 
It is hoped that the results presented here will not only 
be useful to inventory managers and high-level planners 
but will also generate further interest in the application 


of control theory to administrative problems. 
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APPENDIX A 


Z TRANSFORM 


A.l DEFINITIONS AND PROPERTIES 

The solution of difference equations by the z-transform 
Method is very useful, because we can transform difference 
equations into algebraic equations in z. Once solved 
this algebraic equation we can then find the inverse z 
transform of this solution to obtain the solution in the 
time domain. 

mils, Given a discrete function of £, £(KT) 7 the 

z-transform of this function, symbolized by F(z), is defined 


by 


Z[£(kT)] = F(z) & J e(kt)z“ (44) 
k=0 


From this definition we can obtain several useful 
properties. For example, let's find the z-transform of 


£((k+1)T): 


co CO 


Y e((ktl)T)z * = J £(kT)2° 
k=0 k=1 


Z[£(k+1)T) ] ee 


Oo 


= 2[ ) £(kT)z* - £(0)] = 2 F(z) - 2£(0). 
k=0 


aN, 
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Thus when a difference equation is transformed into an 
algebraic equation in z by the z-transform method, the 
initial data are automatically included in the algebraic 
representation. 


Other useful properties are 


zfaf(kT)] = a F(z) (46) 


and 





(47) 


where 1(t) is the unit-step function: 


These two properties permit us to find the z transform 
Mm: a constant, since any constant a can be expressed by 


Bene product al(kT): 
Zlal = az/z-L (48) 
. Another property is the distribution property 
Z[£, (kT) + £ 


(kT)] = zlf, (kT) ]+2lt kT) ] (49) 


2 2 | 


ie 



























The z-transform of the most common functions are given 
in several references listed in the Bibliography [1], [4] 
and [10]. 


The inverse transform of F(z) is £(kT), and is denoted 


Zin ORGZ | = etic) (50) 


From the tables of z-transform we can find the inverse 


of the simplest function. More complicated z-transforms 


may have to be expanded into partial fractions so that the 


tables can be used. Normally, we expand F(z)/z instead of 
F(z) into partial fractions because this leads to functions 
with z in the numerator after we multiply back by z, and the 
functions of z appearing in tables of z transforms usually 
have the factor z in their numerators. 

When the partial fraction expansion does not give 
tabulated functions, we may have to find the appropriate 
z-transform from the definition. This will be the case 
in our problem. 

Similar properties of z-transforms exist for inverse 


z transforms and will be useful: 


s 


a sip(z)/z] = 2 t[2 *F(2)] = £L¢k-1)Ti , (51) 


s gliar(z)] = az “IF(z)] = a£(kT) , (52) 


3Q 





= a4tl(kT) = a howmete = 0 | (53) 


and 


The characteristic equation of a system is equal to 
the denominator of the overall transfer function and the 
Stability of the system can be determined from the location 
of the roots of the characteristic equation [10]. A con- 
dition for stability is that all roots must lie inside the 


Mit circle or 


Another important property of the z transform is the 
Final Value Theorem [10]: 

Wibeex(c) has the 2 transtorm X(z) and X(z) has no 
Meles (roots) outside the unit circle (|Z| < 1, that is the 
Bondition for stability), then the final value of x(t) 


Meor x(k) is given by 
ieee tt) ..= vlam x(k) = Dlimltz—ly xz). (56) 
{70 k oo a>] 


To prove this note that 


Zeek) ) = X(Z), = Bain une 2 
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co 


=I 


Bil scikel) = ZX(z) = Z2x(0) = ) x(k+1) z 
k=0 
Hence 
ex(z) —- z2x(0) - X(z) = (z2-1)X(z) - zx(0) 
= J x(k+l1) ee ) x(k) ze 
k=0 k=0 


from which we obtain 


(z-1)X(z) = 2x(0) + J [x(k+l) - x(k)] z* 


k=0 


Because of the assumed stability condition, we obtain, 


ma z> il, 


which is Equation (56). 
A. 2. Z TRANSFORMS OF EQUATIONS (8) AND (9) 


We found that the reliability immediately before the 


(k+1)th rework is given by equation (57). 
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Since eaves and aR,e 7 are constants, we can apply 


the properties (45), (46), (47), (48) and (49) to find the 


z-transform of both sides: 


=> bt 


AGG) jeez leaner RuRT)| © Ziae ] 


= ZIR((k+1)T)] = (l-ajye ol z{R(kt)}] + aero zc) | 


ZR(z) - ZR(1) = (l-aj)e Ree ae ee 


which iS Equation (12) of Chapter I when we use R(l1) as 
the initial condition. 


Samilaply Lor Bqauation (8): 


Which is Equation (10) of Chapter I. 


m.3. EXPANSION INTO PARTIAL FRACTIONS 


Equations (58) and (60) can be solved for R, (Z) 


Baa R(z), giving: 


QZ 
(61) 
Tee ih ei) Ge bea 
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and 


R(z) a (62) 


z-(l-a)e ie Uiee Wen eT 


In order to have simpler expressions we can proceed now to 
expand into partial fractions. We will exemplify it for the 
last equation only, Since they differ only by a constant 
factor in the last term. 

The first term has already a simple form so that we 
need only to expand the last term. For reasons stated in 


Section A.2, we will expand this term divided by 2z: 


-bT -bT 
iq ae Z = ae 
B [z-(1-a)e >") (2-1) [z- (1-a)e°"] [2-1] 
A B 
= —— + aaa | (62) 
Z-1 Pot oony = Dr 


To determine the coefficient A we can multiply both 


Sides by (z-l1) and evaluate the expression at z = 1: 
en. a B(z~-1) 
=pT se oe! =pT 
z-(l-a)e z=1 z-(l-aje Z=1 


This gives 
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To evaluate B we multiply both sides by [z-(l-a)e 








Ou 


and evaluate at z= Mea ou so that the term with A 
becomes zero: 
wae De 
———. = O+B, 
Z-1 SSG ae lowh 
this gives 
ye PT 
B = t 
bios) eee 
Substituting back into (63) , we have 
-bT -bT 
ae _ ae 
(z-1) [z-(l-a)e Seb BE (l-a)e ie 2 1] 
-bT 
- ae 
-bT 
[ (l=a) =Vi(z—=( =o e ] 


and R(z) can now be put in the form: 


-bT 
R(z) = ———+——-R(1) + << os 
os (lees mee eee yemo 
-bT 
i Zz ae 
eelae\ ones Plo aa noe 


The expansion for Rq (Z) follows immediately from this last 


Bquation: 
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A.4. INVERSE Z TRANSFORM 


The first and last terms of Equations (64) and (65) 
are not found in the tables of z-transform so that we will 
proceed now to develop the expressions of these z-transforms 
from a guessed form of f£(kT). 

Since the first and last terms are identical in terms 
m1 z, differing only by constant factors, we need to do the 
development just for the expression in Zz. 


Thus, let's assume that 


F(z) = » £(kT)z * = » Se ee (66) 
| k=0 k=0 
=e 
Multiplying both sides by (l-a)e z  , we have: 
-bT - S +1 -b(K+ -(k+1 67 
Jee pT Le (2) _ 2 eons 1. b(k 1)T, ( Jie 7) 
=0 
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Subtracting (67) from (66) side by side we get 


k+l] —b (k+1) 7, -(k+1) 


(68) 


(62) 


(70) 


[l= (1-a)e ° 727) F(z) =e) @eo) = O82 = ) (l-a) 
k=0 =O 
or 
Be eae 2 1B (2) =1+ f (1-q) KtL,7b (kt1) TL (1+) 
k=0 
= 3 (L-q) St1,78 (kt+1) TL = (k+1) 
k=0 
or 
HESGeaNGlo es Spica) oa 
Thus 
ale 2 
ee 2 oe 
1-(l-a)e pT. z-(l-a)e bT 
Prom edits tte fOlLlows that 
2 > Zt = E(kT) = (1-a) *e OF? 
z-(l-a)e 
Also, 
| 
zt Z a Gla ea se 
po Nees Gon) Gaal Uae) Cae len 
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and we can get similar expressions by changing the constant 
part. From (52), (54), (69) and (70) we can now find 


the inverse z-transform of equations (64) and (65): 


_ -bT _.\ kK -b(k+1)1 
R(kT) = (1-0) *e 7 R(1) + —2S—_____ + Altra) e 
l-(l-a)e (l-a)e -l 
or 
R(kT) = (1-a)“e *TR(1) + —*2 —1 (1-0) Se PATA 
(l-a)e -1 
ial 
And from (65) 
R_ (kT) = (1-a) *e PTR + ——2 a [ (1-0) te PST} (72) 
(l-a)e Sel 


These are the equations (13) and (11) of Chapter I, respec- 
tively, and constitute the solution of the difference 


equations (57) and (59) in the time domain. 
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APPENDIX B 


COMPUTER PROGRAM DISCUSSION 


In Chapter V we have shown that the functional equa- 
tions of the dynamic programming to compute the optimal 


costs over the set of possible controls are 


CH (R(K) -a(k) ) Je key (ROK) W(kK) ) +98, (R(K+1)) 
and 
Jie (RIK) ) = Se 


The program that performs these computations is 
listed in the Computer Program Section of this thesis and 
is composed of three subroutines and two functions. The 
Subroutine QUANTU computes the discrete values of the state 
Meriable R(k), and of the inputs U(k). For each discrete 
Value, the subroutine STATE computes the value of the state 
Variable at the next state and the dynamic equations of 
Our model (15) and (16) has to be input in this subroutine. 
With the vlaue of R(k+l), the subroutine cost selects the 
Value of the performance measure J by interpolation between 
the two nearest discrete values stored in the array of 
discrete values of optimal cost, COSTOP, of the (k+1)th 


State. This is possible because the computation 1s done 


oy 








backwards from the last stage of operation. The costs of 
the last stage are computed by the Function HN, and the 
subsequent cost for each state and each value of the input 
U are computed by the Function G. These costs are summed by 
the subroutine COST and in the main program the minimum 
value 1s found for all discrete values of U(k). These com- 
putations are repeated for all the discrete values of the 
State variable as can be seen in the flow chart of figure 9. 
The program output presents the optimal costs and corres- 
ponding effectiveness of rework a to go to the final state 
for all discrete values of the reliability so that starting 
with a given initial reliability at the initial state or 
first rework table we can, interpolating through the values 
of the next state until the final state, find out the optimal 
sequence of rework for that particular initial condition. 
The program output for the examples of Chapter V are 


presented in the Computer Output Section of this thesis. 


88 





) is 


Cee oy ee 


eee OM 


eee se IT 


N NX 
20 


s 
+ ” 
— 


tf4 
J) 


JOEMATIVG 


hoe ee 
309 ) 


ed ee PSs 
. iso 


me wa Pe ae eee eee eee ewe eee eee eee ee 


—_ oe SS ee eee eee Ee eee ee ee oe ee eee ee eee 42 oe - 


eer een oe nee oe 


eee ®& @ » © © © © & @~— © & © © gw Oe 
YVerCuiomin wo 1 


Cay POLI VIG: yea ON os UU OG 


eS a ae, Sy a Gy Gp, Ge eg es an, ag a Ga ee a a a 


ay Ud Cis 


en 
Js 


: 


Cer ee 


a 


Oy a Ca ta Oe) Mat) Ge ee, OY oer ay rae ey SD | Ge OE Se fae 
*# ee © ee @ © © » © ®* g@ @ & » * @ 


Gwe OOM OIOCOC OG 


89 


Le ee) ate 





Sey Tee oT Cut GxXse 


PA Saas Vo Ged eB) << 


Bun. ENT NixrT Ses 4 Sea Ok. 
ois, | Beant cat nl, oS! REWORK 


== eee ca cabo <a 228 cm oo ce cee eee eee ee ee ee ee eee eee eee ee eee ee eee ee ee 


1.U0 ATS lors Je% ve Si) 


=e ee eee SS RP cee we De ee eee ee ee eee ee eo ee ee ei ee eee eee ee eee eee oe le 


iano JIe7D4 e+ Car 


het ee Peat 2 S22 ee es SE ewe OP eee ee eee eee oe ee Oe See ee ee eee ee eee eee ee ee ee ee ee eee en ee ee oe 


eas Je 247 Der Orie re 


Ue St Ue 3s Jer CPase 


eat SSeS eee ee eee ee eee eee eee eee ee eee ee ee eee eee eee ee ee ee 


ere Geos 2 ‘J ee Orcs, 


= oe a ee ee ee ee ee ee ees a ee ee eee ee ee ee ee a eee eee ee ee ee eee ee eee a 


nm “a 4 5 “ ¢ fy 2 
Uonenedt oe: ee Wet w 5) 


=e ae ee Oe ee ee ee ee eee eee ee ee ee ee eee eee ee eee eee ee ee ee ee eee eee eee ee ee 


Deer e523 e+ Greco 


o_o 2D ee eee oe ee ee eee ee ee eee oe oe eee eee ee ee eee eee oe 


‘ SSNS Ries \ ee ‘ = rr ¢ & 
<de. 2) os J es WI eae wedi 


= ee ee ee ee ee ee ee ees ee ee ee eee ee ee eee ee ee eee ee ee eee oe ee ee ee eee ee 2 


es nS Pe Oe S Gece 


oe CP ee ee eee ee eee ee eee eee eee ee ee ee ee eee ee eee ee ee i ee ee ee eee ee ee ee 


Oey el de> UPA 


i Cie eo ee Gacy 


—= SP wee eS) OP OE es S82 ae ee Fe eee ee SS eee eee eee eee ee ee oe ee ee eee eee eee ee Se oe 


7 VAD ry fa , 
ea Jet? e* Geee 


J .4 0.40 


pa) 
_en ween SP ate Ses Eee oD ae eee oe ee oe ee ee ee ee 


ves 2 Wisma Je + Jeo o 


Aen, | -m 4 F, be he I NX i. cs 
EC my 0 ee Ue rt Greta 


= aD ee ee ee eee es ee eee ee ee ee ee eee ee ee eee ee ee eee Se ee Oe eae ee oD cee ~<a 


oe yess U0 Jie + cen, 


— oe oe ee ee ee ee Se ee ee ee eee ee eee ee eee es ee ee ee es es a ee eee ee ee eee ee ee eee es SS ee 2 oe 


90 





eee oe er a) ee 7 Case si 


—- 


reise ale te eee immed Cs 2 atl ests ANI) cor 


-— <e Se Se ee oe Ge ee ee a ee SS ee ee ee oP ee Oe ee ae Se ee ee eee eee eee ee eee ee oe 


Paseo Tall fle Xx), MINT “UM Tey Py ai 
mo( kK) Ao(K+te) Cars RENGRK 


ane ee ee eee ee Se ee ee ee ee ee eee ee eee SS ee ee ee ee ee eee ee ee ee ee ee oe eee ee ee 


£-CC Oe se Cat 0-90 


—_— ee oe Se ee ee ee ee ee ee eee ee ee ee ee ee eee eee ee ee oe 8 


Ore 73 0. 753 On Uc 


ee ewes = 2p ce ee ee Se ae ee oe Se eee ee ee ee eee ee ee ee Se oe SP 2 eee 


Ce oo eG eG Cae Cao? 


—_ a SP a ee eee ee eee ee ee eee ee ee ee ee er eee eee Se eee ee ee ee 


ae <Se ae See ee eee Se ee eee eee ee ee ee ae ee eee eee ee ee Ee eee oe eee oe 


Cen Sa oe 747 Sie t Heh 


Oe ee ee ee ee ee ed ee ee 2 ee ee ee ee ee ee a Se ee - 


s 


Ue ss QO. 743 oan recat 


o . ~~ A ~ 7 {fy a 
~ + - — “~ 
Cielo + OC ea eee Sad Oe Wwe Sk 
2m SP ee Se ee eee ae ep ep PU. le UUme USP le Er Cee US USE Sr er Oe See eee eweamw oe «ap am “HB ece SS 2S = Be oF «aes 2 


— SP ee CE Se eee a eee eee See we eee eee ee ee ee ee eee ee ee ee eee ee eee ee ee eee ee Se eee Se eee le 


Or te epee Ce? Je 5Y 


ot 
—_— a ee eee ee ee ee 2S eee ee ee ab a oe SP Se ee ee ee ee ee ee eee 


a = . ~ mY an , 
ee | eon 3 ears ose 


ere GE 22 2 4 SES 2S ees aes SES See SP SSE EE aR Sew eee See ee SS eee ee eee OS eee SSE Se eer ae ee ome ew ee Oe 


bec Je 444 DEA 6.5) 


=e Se eee oe eee ce ee ee ee ee ee ee ee ee ee eee ee ee eee ee eer eee eee ee Ser Te! Se Se ee Se el 


Je fs OG. 920 Gian! Jie Su 


ee | Ge Flo Cat Geog 


— =e em Ee ee ce ee ee es eee sie ee ee es ee i eee ee ee ee ee es ee ee ee 2 ee eee ee ee ee 


os as OF a Cay wie 


caves Ge a Oef U.0V 


Ge 44 Jems Ge f Sie OL 


—_ = am «8 ee ee eee ee Oe Se ae ee ee ee OS a ee Ge SE SE ae Ee ee oe Ee oe eae oe Oe ae SP ow ae eee eee SO 


OR aes Oe 2 Us Digit Bere, 


ee cee ee eee oe ce Se ee ee ee ee ee ee oe ee ee eee eee ee Ee ee eee EE ee Se 2 aa es oe = =e 


ee 3: oe eas Giant We a 


—_ = cee Gee oe ee ee ee es es ees ee eee eee eee ee eee ee eee eee ee ee ee ee 


fs 


Gere! Ue oe Je ? 0-60 


=e ee ee ee ee ee eee ee ee aa oe ee eae em ee Sa eee Ee Jee ae SP oe Oe ee ee ee ee eee 


Os 8, Di a She Qe 0 


ope 





Cie Jee ih oul Gy eue 


Sirs ae Oe ere GC. 50 


=e 2 eee 2 «se ee eee eee See ee ee ere eee eee ee ee eee ee eo ee eee ee eee eo 


i Je 75 ioe ue EY 


dP) Oe ou rent Cs) 


Wie oe Se a. ions ie 0 


Bly a Gage 25 et ate 3) 


=e 2 eee ee ee ae ee ee See ee Se Se eee eee eee eae eee eee ee es ee eee ee ee ee ee 


2 
. 

—_ *t 
diet.) 


—_ =p ea 22 ee eee ee Ee ee a a oe 


re | 
Cu 
6 
i) 
ti? 


ar} 4 [~ , 
Py Lue es Se rere, 


eb EE Pe a6 ee 2 2 awe es See ee Se ee ee Se eee ee ee 


— a a 2 ee ee ep eee oe ee ee ee ee eee ee ee ee ee ee eee ee ee eee es ee ee eee eee ee eee ee = 


Jeo Cs oo Lae Ce ow 


Cea | aos ‘ee arc 


> oe ee ewe eee ee ee ee ee ee ee ee eee es eee eee ee ee eee eee ee ee ee ee eee ee eee ee oe 


Seals ihe oS be lect Ceo 


0 tore, Oe FU dea Ue cw 


Je S54 On205 Lad Ce, 


me CS ee ee ee ee oe ee oe oe ee ee ee ee ee oe ee oe eee eee eee ee ee ee ee ee ee ee ee ee oe 


F _ i 4 cf 
Deo S Pn Oa Lae Wa a 


ewe wea ees ae ee ee oe Oe eee oe ee Se ese Ee eee ee eee ee eee Se Se KP eee ee ee eee eee ee oe 


Je52 0D. ote / Led GC. 2V 


—_m oe ee SEC ewe se ee ee ee eee ee ee ee eee ee Se eae eee eee oe ee ee eee ee ee ee ee ee ee ee 


ope Dec Gen Ue 30 


JeEU ee ow ans G.60 


we a Se ae oe ee 2 ee Soe ee 2S eee SS Oe es es ee ee Se ee ee ee ee eee eee eee ee eee 


ge 





ae eee Ce geil yee ish 


ee N fF he re Tr A PT 
Poth) Aste gost aS 


1201 Die oe 2 eae OP AeeTS, 


Ce 4 G23 54 Ve + JVerJ) 


—_ om Se eee ee Pee eee eS ee ee eee ee ee ee ee eee eee eee eee oe eee ee eee ee ee ee 


3 = Pa 9 p ‘ » 
CaS Oe” Ce eta vecu 


—_e eee SE ew ee 2 eee we eee Ee See ee ee See ee ee mcr eel eee eee ee eee eee ee eee eee ee Se Sele eee ee eee 


Cees Le 247 eee gt) 


—_— =P ae a ee Oe ee Se ee ee eee eee eee eee ee eee eee ee ee ee eee ee 


7+ Ue: 


= TS eee eee eee ee Se eee eee eee lee Seer eee eee eee ee ee eee eee ee ee ee eee ee eee eee 


eo Oe 9 23 ce. Us 31) 


oe A es ee eee eee eee eee eee ee ee eee eee eee eee ee eee eee eee eee ee ee 


bg Ss: Js Veet Ola 0 


=e ee ee eee 8 ee eee ee eee eee ee eee ee eee eee ee eee ee eee ee ee 


en cue Heese We Oe oU 


=p a SP ee ee ee ee ee 


Ceo Oa 1.4 Oe, 


—_— 2 eee ee SSeS ee See eee Ee ee eee eee eee lee lee ee eee oe ee ee eee ee eae ae ee 


ence. Ge 20 tert veo 


ee ao Ono Vert Oe 30 


ae ep eee ee oe ee ee eee ee ee eee eee eee ee ee eee ee eee eee eee eee eee 


Geos Je IU? lo. + 0.259 


me oe ee eee ee eee 2 2 2 SS ee ste ee ee ce ee ee ee eee eee eee eee ee ee ee Se 


Ceo gS 05 Let Serre 


ap oe ee ee ee eee ee ee ee ee ee eee ee eee ee ee eee ee ee 


Pn i! = ANY : Ze 
2 Pate Wren ieee OPS), 


=e ee ee ee ee ee a eee ee eee ee ee ee ee eee ee 2 ee eee ee oe oe Se eee oe 


Gen oa oo) a DoS’) 


a = en eee eee ee ee ee ee es es ee ee ee eee ee ee eee ee ee ee ee ee oe 


Crane Je ras 1.4 O.5U 


‘ 
= 
Fe eT eT i i i i i 


Cees 0.3 70 4 os Oe 


ee ee eee ee ee eee ee ee aes ee ce ee eee eee ee ee ee eee ee ee oe oe 


ce 





Ol 


SOG GAT oe 


Se Er 


OA TA INPUT 


le Seca rt on 


ce J La 


~~ 
' 
—_ 


TPA Ay ac 


ava, 


fe See 


FIN 


woe be ade 


~ 
{ 


ree ae er 


“TNT AU" 


DO aa te ay 
Ns > im Ye > | Tew LP cm an) Bs A) 1 Pon 1 ah as | 
oe ef @ ee ©» & © eo © & © & © wo © F @ 


DO OE HAF AOL ON OID IE AIEO 
DOr Dal Gi cy CC Cie iy Ol Carr Ce 1 


me IOC OD) DO) OD OO ae Oy 


94 





Satya eR SU hey Se SHS 


TASLeE FOR THE 47TH REWIRK 


Pe meer NEXT MEO “Ui OPT AVAL 
erat K} RS(K+2) Gast RE WORK 


= 22 @m ao) ce om cee cep cme es em cer eee SR om eee Geemee cer we om cer cer ae wee Sm ow SP oe oe ee ee eer eee oe ce 


t ro e939 Osu Gece 


Rt es 
CeS5 O-e990 Gren Je 50 
me me ee ee eee ee ee ee ce ee ee ee es es ee ee ee ee ee ee ee oe 
oe a ee ee ee ee ee ee se ee ee oe ee oe er ee ee es Se es ge ew we ce 1 ~< 
I TT eee es ee 
De wes Je G7 er 2 One 8) 
ee em os ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee 8 = 
en ome re ee ee ee ce ee ee ee a ee ee a ee a ee oe 
—e me ee ee ee ee we we ee ee ee wn ee ew ee re ee ew a ee a ee 
meee ee ee ee ee ee ee ewe we oe ew ee ee es oe es = 
See ee ee ee ee | 
— cee a ee eee eee ee ee ee ee ee ee ee ee ee ee ee es ee ee ee aa = 
eee cme SE eee ee ee ee ee ee eee ee ee ee ee ee oe a 


—_—m= “SB 2S eeeeeanwe Fiweenw sew w ease eae & 2 Fan 68 Sse SS See a ae wwe eS] 2 as Gs = man ate See eee ee eee ee 


an =) GE SE aS Sea «st a «ee ome See eee Gee ee een ee eee ee Oe a ee ee ee ee eee eee eee Se 


O.35 DOE ee are een) GeP) 
G.°34 bev lf Cres, Dg) 
Weds eo tS ec) Jae 


Ee SP SE aewenw amas SP eo Eee = GE Se eae Se eee Eee ee eee ee See er ee ee ee ee eee eee eee ee ee ee ee oe 


J6 8: Jeo 74 Jed Jes 


—— pe ee ee eee eee ee eee eee eee eee ee Se eee eee eee eee eee eee ee ere eee ee ewe oe 


== SP aes ee Se ee es ee ees es eee eee ee eee oe eo ee ee ee ee ee ee ee ee ieee ee 


25 





Coe ITER ATU GASE 2 


Srl NOS MINT 4U 4 NP AMAL 
K ) 285 0K + ]) CSS § 2E WORK 


ae) a) UeISU Oe J Ue SU 


—e——— SS SS a ey oe OD ee se ee ee ee es es age Sas mam aD = ames == 8 oe ee 


Se CM em eee = a “See eee SS SS ee om cee eo oe ea “Sob SS eee eee SS ee cee Sew =e oy EP eo Wee “Se oa =s eee = se op = 


== a cee oe oe oo — ew cee ee ee ee ee = ae ee ee es ec =— ae a SS SS oe ee oe a ee et ee ee ee 


Creel OPuatr eG oe Cao 


2 oe oe oe we eee Se Oe ee ae =e — ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ae ee ae Seem 2 a= 
ee ee ee oe ee SP ae ae ee ee ee =— et ae ae Coe a —— SE CRP ee ee ee eee ee ep oe eee ee ee oe 


oe ce eee ee «< ee ee ee | ce SS ee Ce ea cee ee | = as oe —— a ee oe Oe ee ee eee oe oe oe oe om a 


96 





TOY4PITER JUTPUT GNSEs 


sea) EAT MIA TMU 4 Nie tts 
(K) Stings 12) (SoSe8 REWIRK 


—_a a a a ae ae ce a ee ee eee ee eee ee eee eee ee ee eee ee eee ee eee ee eee eee 


=— om Se SP ee eee ee eee ee ee ee ee ee ee ee ee ee ee ee es eee ee ee ee ee ee ee ee ee ee ee ee ee 


On Fe Ce 7 1 Qed Ue 5U 


—_ oa GP eee eee ee ee ee es ee es es ee ee ee es es eee ee ee ee ee ee eee ee ee ee ee ae 
> eb see oe eee ce ee ee eee ee ee eo ee es ee eee ee es eee ee ee ee eee ee ee ee ee ee ee ee 
= am ae ee Se ee SS ee Se ee eee oe ee ee es es ee eee ee es ee see ee ee se ee eee ee eo 
qneanm eat aut SE we a ae ae aw SE C2 ae ae oe? SS Aeaweanm RP ame eee ee eee ee eee ee eee eee eee ee 
qn eam Gece aw ep aw 42 SS eee 2am 22 eae aw Se CP GP SE 2 Ee Se 2.2 SE ew aE ee eS oe ee ee ee ee ee ee 


— ams 2 Wee Se ee ee ee ee ee Sa ee ee ee eee ee ee a eee ee ee ee ee ee ee es eee ee eo 


Oe 72 Jet 35 Jed Coad 


amt Gpeeam a Se ee ee eee ee ee SE ae ee ae ee ae we OOS SE ae ean OR Ce oO Se SE OE Se ae eae mim ee eS DP 


reac 2 aL 9.0 Ge SJ 


ee ae Ge ae Gee ee ee ee es ee ee ee See eee ee ae ee eee See eee eee a eee eee ee ee eee 


Gers 2 Je 98U Ue J OP ese) 


a oe ee Se ee ee ee ee ee ee es ee ee eee eee ee ee ee oe ee ese ee eee ee ee Oe eee eee ee 


a5 f Cie 1s Jed et) 


ee eee eee ee ee ee eee ee eee eee a ee ee Os Ss oe eet Bee Ne SP mem eee ee eee ee ee 


os co ee ees ee ee ee ee ee ee ee ee ee ee ee ee eee ee eee es ee eee eee ee ee es eee ee ee ee 


= ee ee es ees ee ee ee ee ee ee ee ee ee nee es ee eee eee eee eee qm US ea! a am ae oe eae we CER oe Se Ss Se 


O06 > U.S 79 Jed GetJ 


= ee eee Swe se oe Se ee ee ee eee ese ee eee ee eee eee ee eee Se ee eee eee eee ee 


Cece 9.974 Jed ee ree 


a ee ee es ees ee ee ee ee ee ee ee ee ee ee ee eee se eee ee eee ee ee eee 
eee ee ee es es ee ee oe ee ee ee ee ee ee ee ee es eee eee eee ee ee eee eee eee es ee 


Ue Eu Gea te QV Ca 2) 


ee ees ee ee ee ee ee ee ee ee ee ee ee es ee ee sa ee ee ee ee eee es eee eee ee 


OF, 





Ce ee eee Uy Case 2 


Peomeeer CF Hee. |r REnWoaKk 


Ee apam “PSP am awa EP am 8242 ae aw 64 Geet ae SR ae =D aap a 4 ans amp EE GE 4 we oe a Ge aap at a= atic .27 a = a 


— = ae a Se ee ee a ee ee eee ee oe eee ae a eee ee eee eS ee ee ee ee ee ee eee Oe ee 


—— = ae See eee ee SS ee ee ne ee ee ee se eee eee eee eo ee ee ee ee ee 


Ear amp ae EP (Ri ea Ee “EP eee SS ae awa Se SS See 4 SEO AE SE OEE Se SS SS EE SE ee ae ae 2 ae 2 ae a ee aa 


eee a eee eee ee ae ee eee eee eee em aoe se eee ae ae eee ee ee eee a ee ee eee ee ee ee ee a 


Ep awe eae =e eae ewase oawemwe= Baw ap ss ae SSFP ee = oP ae See ee SS ee ee ae ee ee ee ee ee ae eee oO eS a oe 


== a a 6 Ge 6 SS ee ee ee eae ee ee eee ee es eee eee ee ee ee ee ee eee eee ae Se ee ee oe 


Deas oe 1357 ue) Care 8 


i «at “Peesease ae ii aeeep aeieanam 2 eenBaewaweam ae at 28 ae baw ae ae amwamam 4A Ge ae -S ata iat aan 42 282 a 


=> Ese Gees <a ae See cee Ge ee ee ee eee ee eet Set iO ee eee eee ee ee ee ee eee ee ee ee eee ee 


=—— Gp Gp Cr ae ae ap wee ae ee oe ee eee eee ee ieee eee ee eee ee ae Eee ee ie eee ee eee ee eee 


= aap «ap am GP see eee == cep i ee cee oe se ce Go see oe ee es ee eee ee ee ee ee ee eee ee ieee ee Oe eee eee 


beac 02S 74 Ged Ge 30 


= 2 a Oe oe oe aD Se ee eee es eee Se oe ee ee eo ee ee es a a ee ae ee ee ee es ee eee es ee oe ee ee ee 





OAC Wee Sue UT 


7 


OA tA Tee 


xc) 
“{[ e 
TFT sf 
~ 

2 
— @ 
oe 
< 


aie 


~~ — 
ot 


J? 7 tere 


1o 


A 
“4, 


ats 


ap SOS ean, SF] eww See wp ee SS aaa es eae We a ee Se 


—_—EPEae SP ape Pe ae Se ee ee eee ee eee ee eee Se 


IND OF9N4 89 AHO DED OW Ptah Ht4 
DMAODCUOT POPPIN OPH O Hb Abe 
oeee ff @ €& @ @¢ we © & © & @© @ e@ 8 @ 
Na Fee at Ot a is co Gay Int om 
Veet 4 


DO IE HO SR be TN OD 2H 8D 
Th DONWON Et ODO MN MMM wen m 
e @eeee#e#*#ee#ee#e#e #*€ @¢ #¢ @ @ 


(DP can am OD am ap lee Pap am an lamas alee 


— ae ow SP cee cD seb ee eee SSP Se eee ee ee eee eee ee 


oe) 





CY4PUTES JTPUT 


Cy 
a? 


is) 
Mm 
Wo) 


ee roo eS he SHREW AK 


COE en NEA ENA JPT MAL 
mo Cho) NS ne iD, COs seal eins 


—— ae 4D GER om ae ote ap a op Gem ee EP Ree ee ee eee ee ce ee“ Pe ee eee ee ce oe le 


1.UG Vevisz QO. 5 cai 


—_ SES amp GED GS a om ae aa 188 EP oe a oe oe WE <P ae oe Ea ee ae Wa EE oe oe eee oR es oe ae fe eS Se 


Cero: 02999 C. 5 pees 


—— = GR at oe «<2 M8 Hew SS es See ee Eee Oe eee Se ee ee eee eee ee eee ee oe 


— ap «DP op ome oe 26 ee ewe ee Se ee ee ee eee ee eee ee ee ee ee ee ee ee ee et et 
—, ™ ~ os pw) —_— 51 Sa Te 
{ Ma 
Ue - | eof: el Jed) 
ae aes se awaeamanmweamwranweeawawenmwran ew ae aE SSB ee See ESS SS SS eee ate ae ae EE SP ae a ae «me Se SE ome GP aw 
. ayes \ Sb) = ve ; . 
Wwe 4% we is Ue (ee + 
=e ee ae ee ee ee eee ee eee eee ee eee eo ee eee ee ee Se Se eee eee ee eee 
a =~ = I a —_ s SINT cet 
— — : 
ers Ue 73 wef \/ © > 2 
= oe ee ae ee ew ee ome ee ee ee eee ee eee eee ee ee eee eee ee ee eee 


mem <P - Gee a eee ee ee eee eee ee ee eee eee eee eee eee ee ee eee 
_ * es Ian) r ~y ~ =~ 
“4 - ; 
Je 5 Oe i Ome, Ly 6 4 

o—e oe ce 8 eee oe ae ee SP ee ee ee ee es ee ee ee ee eee ee ees ee ee eee eee es eee ee ee es ee a ee 


io | 
id @ at 


Jer Ors Sie 3 5 


oe eee ee ee ees ee ee eee ee eee ee ee eee ee ee ee ee ee oe 


=—_ a ob ae a 


oA Ge 23) S28 Gane 


= ae aD Se ee Oe eee ee ee ee ee ee es eee ee ee ee eee eee ee ee ee ee 


‘ So lk ae la Be id ta aN ~ _-™; 
Wei lt ae a Se Ons 6. 

= Ge 4 ape ae SS 2 Se SP ee ee ee ee es ee see eee oe =e = a= eee a ae <P eee cee ee cee 4 Ge ee SE ee oe —_— oa ae ass 
J ed * Ue ONE =) Le 


ee oP eee eee ee ee ee ee ee eee ee 


hwo 
ion 
ic 


Jai OQ.2°79 Oe Ve Fi 


“amp <9 °38 eer ae Se es a ee eee ee eee eee ee eee eee eee ee oe eee eee eee eee eee eee a ee 


OF Rea) Oeics oe Ve FU 


= , a7? yoo yO 
ue 3D We tai vers ee. 2 


ee -  ee  e e e eei e ee e eee e eee e e e ee eee eee e 


0.354 Olt fe OLS Oe ere 


==> > a S65 ap Ge cD GD come G8 A ee ee ee ee eee eee eee oe ee ee ee ee ee ee ee ee ee ee ee 


nes Ueo t> 0.9 Cas 


= aD cop ac «aD map Pp apap asamp (ean sS SSP cepa <P eee eee eee SP 6 cee ee ee eee ee eee eee ee eee ee ee eee oe 


wy eee 92974 Ceo Ceo 


[> aD ce Pee GED SE ee a ee ee Ce ee ee ee eee eee eee ee sees eee eee ee eee ee 


ge Gl Oy eee 0.9 Ca 


<> <GEe <e  ee  e e  e e eee eee  e e e e e eeeeee eee ee ee e eee 


Sah, Dien 2 Cee Jase 


EE SS AS ce ee ee ee eee ee ee ee oe 


100 





COMPUTER JUTP IT CLS re ae, 


(eet h sa Pima Shr 2 Samii 


Cun tc NT “eX MiP i sh jo 
= SO 2S(nr+!}) GorRSa) 2 


a@pa—wt-«2 46 G2 S82 sae ae® 2p ese ae == SP ot Oe ese aie ae OP ast at SE ap ae SE = =P) ae ee eee SP ee eee ee eee OS P aw 


eC Ges e2 ee) Oe 2D) 


eapneanbas =P 442 2832 22 -e ew e8 aoe SS ap = ae =e ee = 2 ee eee SS ee eee SS Oe we ee ee Ee eee OS EPP eae ase = == = = «= 


Ceo ie Peres ee eou 


eq_2E2epanr 2 APap Ss aan Se 4eaeane eae aw w7 at Se se sweaeese oe aia Fe ae = SP ae ae «ee awas 4-8 ee 6 oe SR oe ee OS le 
a aon “~ ito a 
ae 8 We F742 eo) Deco 
qeoeab ase eee an Sant aa ae 42 42 42 .s qpem 42>-op @@pew 22 a2. 26a == 2 at 22 2 = == = pYQE Ge ae FSR ae 2s ae as == sa 


Q.o4 OATES & ea ee 


=—__ ape OP 8 28 ae Bae ae Ee Eee ee eee eee eee ee eS eee Pee el ee aenp == .2 82 == =» = —-_ 


Ge oi a ss ae ens 


eat eae eae oS Pe Tee eee See eee eee ee eee Oe SP aee SE aee aee eeeee eo a ee ee St SP ee ae e2 e e e  e 


rae See, 


—_ 
eee 2 PP Fas SFastbast SBS we EP QF ata BF ab Sas ae aa SS = 


Ved! Vet ss re Pee 


aap abat@esp at 2 aE ae ae aes aes SP eee SP 2B ae MR ata as See ae Pee Oe ee ae ae ese SEF SS ee oe Se Se Te 


ae 2 ie a as 7 . a 
ic ; = 
we 37 Je Sinn te Jc eget s 
amen QE 4am ant EP ane 2a 4 2 aha 42 «s = SE eee ae aD oo Eanes 2Re ewes 8 ase aeiae ae Se ae 8 ase awa SE -—=s Be ap = =a 


es) 9.2315 gel De 


"> 3 t erie I > _ . 
eve eg eae Neier Geo 


Seep ean Baw ae we ape Seat aS ae ae Ce 2 ee Ree a eS eae eee ee ee eep eS e SPee E ae 6S Se ae Se SS a 


Jens ier eet aos 
Seno Jeet ae ee 


re: 2 Die 2 or esate eo 


o_o ae ee oe eee ee ee ee eee eee ene eee eee oe oe wee ee ee ee Oe oe eae S22 =p 


We aie aS ; Oy ae 
MJ ee Voge vee! Oe 


=e 2 =) 2 es SAP et ce Re asa ae SP Poe SP ee ae ase awa ast =] oD 


Ie 2250 tei Jedd 


lat 4S ates 26s ae ae SB eweti es 22 oe 2 ess Se Pe ee ee a ee a oe ee ee eee ee oe Ss See =e oe 22 ow 


LOL 





CHePoTe <a reT ONS 2: 


TAGLE FUR THE 2TH 2E4IRK 


ee nec YINI MUM JPTIVAL 
RS(K) ee) ana: 2EWURK 


CO I a Rt ee ee 


1.00 Wie ee bet Je 30 


meee See ee we FS eeu 22 een ap 82 sea =e EE shee oe ee ee Ee Ses OU a Ses. oe 8 ee STS See Um mevlUmCe lO eee eS ee eel 


OR ee, 32953 le er ed 


_eEPu “Se =e Bee eee 4S ee SS eee se ee eee eee eee Ke eee See EE PL Pee ee Se ae ee ee ee 


Ce Se Uren 1.4 oppere, 


=—e PS Ee ae See OEP ee ae eee ee eee EE Ee ee ee ee ee Se. See See oe SF eee Ue Pl Ce Um mele Se le ss > —.» = .» 


Geol Oe 250 Le Carer), 


Se eee, ee ee ee ee ee 


Ce c¢ Ce eter Ngee ora 


= oP oe eee eee ee ee eee SE ee eee eee See eee ee eee Pee ee eee SS Ue SP lee le eee 


Ce 54 e739 1.+ OEp ye, 


ooo Se le + J.5) 


ee eS SS Pe ee 2 Oe OR amu en Gee fe Se Ee See oe ee eS UD 


SO Ste Gea Ve. + Je 3) 


—— oP ae See ee ee ee Se ee ee eee eee eee ee Se eae ee Se Se eee eee ee ee ee 


Ceol Uae os Oe350 


=~ ee ee ee ee ae ee eee eo eee eee eee ee ee eee ere eee eee 


Cec < 0.92% Le + Ged) 


ei =e P22 ase ewe we ee Ean =2ae = == aa «- aeweegp any 83 F' 2 Ss SCP wae ee ewe ae Se =F 3 SP ae aerw 24 ane Ss = oe == 


4 


ieee Gar 2u tet On) 


— ean ao 28S aS OEE OE ee awe ae OP QS eae 22 ae 26 SP =*s SBS Sea ae ae See ee Ue SE ee 6 em awe a Ga «a 


Gragert Osols Le + Sie 


—_ we eS eee Se ea ee ae ae OP Se Oe oe 2S eee ee ee eee ee ee ee eee eee ee ee 


Ge 85 Oe Ler Ja50 


Oaks, oe oc. 1. + JS 


Dans Oe. 70: 1. 0 ecJ 


ue 82 errs Let Use 


ae oe ee ee ee eee ee eee eee eee ee ee ee ee eee ee ee eee ee eee ee ee ee 


C. 3h Wignoas Let Je dJ 


ao ee oe eae ee eee ee ee eo ee ee eee eee ee eee eee eee ee 


Oren! eo Te + O22 


ee Sea SSP ee Eee ae ae es ae ee SE Oe Bae ee eee See ee ae ae eee Se SP. ee 2 ee eae ae ae SS a 


HO 





Soe OTe a. ree Coe. «2 


CURLZENT Wer MENT MU’ cone 
eG) mS hy Cre 5p : 


—= ame EP WS D> ee ee eee ee ee ee oe ee eee eee ie eee ee _eomeaw 4 Ee ewe ae ee ee awe 4 42 ae as as GE 


ION) Spee a eet G.e33 


ae 22 an 22 eee =e ewe See SS eee owe ee eee ee ee ce Se Se ee ee oe ee a ee ee eee ee ee ee ee ee es ee ee ees 


Ge 35 Je 32 oe Le 50 


oes 4 Gis = aes ef ee a.) 


= eee eS eee SD 2 ee ee ese eee ee ee ee eee eee ee ee eee ee eee ee ee ee ee ee ee ee oe = =m «as 


ry ~— a ~ so + . Seay 
EA er oP Ve a 2) wee, ee 


aoe ee eee ee OE Se a ee ee ee ee eee ee es ee eee ee ee es eee ee eae ee ee ee eee Se ee ee 


‘ ~ 4 . 4 »* a 
ex 4 OSs eS, Gare 


ene ( @2eawee= se awe eae 2a ae 282 ae ae 2 aS 22am ew eee ee ea eee eS ee ee eee ee ee ce ee ee ee ee ee ee eee ee 


’ F in me Dal “4 a hia r- ? 2 ~ 
i= 3 Wes oo va Le OU 


eo Ue oe. Pays Gaeo 


_ ae 2 ae ee ee ee eee SE Se See ee Seem ee ae Oe eo eS! eee eee ee ee eee ee ee Se 


3 feet net ees saat t Gale 


[eee «ao Soe 2 ee 4 2S 2S Se eee ese Se OS SE eee ee Se ee ee es es ee eee ee ee eee ee 


Oye ee oe ae ier Jee 


a 2 ee eee ee ee eee 2 Se eee Se 28S Ss 2s 22 SS Ss SS Eee eee Te SS Se So See 2 ee ee ee ee ee =e 


Je 5c Ce 320 end, Ge J 


= oe Se ee ee ee ee Se eee ee eee ee ee eee coe eo 8 ee os eer eee ee ieee Se eee ee SS ee a 


os Sot ie AG in 7 Ceo 


new Ge 22 e= 458 2 "3 a eo a 2 aes 2 OO See Cee aw ane SS SS 2 eee Se ee ee SS Pee eee 2 Se eee ee ee 


; a ad ’ ay ct vs Fe —~ - 
Oe oc SS gi een ie Wwe Ti 


Damp 22a 2 eae ea ae ee 2 ae 2 Oe ee ee Sp Se Ce oe ee ee oe st ee eee eel eee 


Cares Qe2ue vet emo 


we ap S22 ae 2 ae! eee -8 22 eee =—8 242 42 oS SF eee eS eee TS Se Ee ee SE oo Seon On ae Po ee eee 


3234 ers ae Oreo 


=e aw ee ae a 8) Oe we on ae oe eee Saw ee SOE ee Oe oe ee ee 


eke 2 Ue Ful ted Jed) 


—_—— a oe OE eee ee es ee ee ee ee ee ee eee ee ee ee ee ee ee ee Se ae ee oe oe eee eee 


Jed 2 Jes ai Ln vec 


a ee se ee ee ee oe See ee 2) Se ee ee eee eC eee eee ee eo ee Se ae Sets = 3 oe Se ee SS 


Sipe! Miers aU et Ure, 


Or 





COMPUTER UT PUT LASE 4 


UAT& INPUT 


‘I MX NU MOGE 
“ 20 30 U 


AMIN XMAX UMWIN tua 
OeiG ewe Us omens 


JF VOSR ATION 4819 ST ue ery ale aye. 


ee erent Cees Tae 


~4{ 
i- 

eo) 
r- 
(Ti 
a | 
ey 


FI aL *INIMG 
Saye Siog 
be 00 ves 
Jey? eco, 
UY. 2.06 
Qe =7 Aa 
Dos Oey sy 
err) Qe 21 
J. 34 gree 
0.93 = S08 
Je CF mS 
Qe Fk Sia, 
eon ge & oD 
eo G ar on, 
Oe sft Sei 5 
Wend 0 
Oe a> PEA Gis 
0 CS ces L+.?7 
eed Tain 
Jest) eS 


104 





CJAPUTER DITPUT CASE 4 


ue ia SI Mel= 2) apes (Olas 


CUTRENT NEAT MIMI Mts‘ OPTIMAL 
RS(K) 25 (K+) Bas 5 REWORK 


anew ee GE ee ae ee eee ee ee eee ee eee Ce ee Cee eee eee ee Ce eee eee ee eee Oe SUE ee eee 2) ee a ae oe 


t 


22 OC J. 9cY Jaw Ge 3S 


=e ee SP TE ee eee oe ee ee ee ee ees ee eemwant GE amwcre Se ew aw =F 26 ee ane ew GE GP SP i. Pema Ge =e ae 


owl O23 oi be 30 


a a SP ge oe ee eee ee eee eee ee ee ee ee ae ee ee eee ee eel eee Pe et ee oo 


Ca 2 Ue 7437 Jedi Jace 


ear Pawar, 2 ae ea= a= Ge aw 2 2 See CE SEE oe eee eee Se ee 2a Oe ee Ge ioe Fe oe ep 2a ap 42 «> Ge 


oe ome 1. < 
wl @ ST we TO Oise Le 


cx = A ae ad i an -~ 
Ce SH ite too aur Ge eu 


-~ oF s = x 
SR ialas 2 ea Ge DU 


[Gor 

e 
{oO 
isd 

Cy 


~ 7 Veen oon rom = , { ae 
Ceo? GOe3953 Cue aes, 


— = oe ee eee Se ee ee ee eee eee ee ae eee oe ie eee ee ee eee ee oe eee ee 


Vig eae Ue gel Jet Cie 


ee ee oe ee eae ee ee eee ee ee ee ee eee ee eee eee ee eee Oe Ee Se ee ee ae ee 
. 


hy ae be Oa ee re: ees 6. 


a a ee ee eee eee ee eee ee ee eee ee ee eee eee ere eee ee eee ee eee eee 


Jeo ft ve Fis J. 2 ie 


=e cee Gee PG Se oe cee cee ee ee ee eee ee ee ee ee ee eee ee eee ee eee ee eee eee eee ee 


oes Oy 3S ae: Cece 


—- ob Go Soe oo OP eee Oe ee ee ee eee ee eee ee ee eee ee ee ee ee ee ee ee eee eee ee ee 


Ue 55 Oem hee Ge ou 


Poa 2 eee ee ee ee eee eee ee ee ee ees ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee eee ee 


Sera st aod ee Cae 


re a ee eee eee oe Oe Ce ee ee ee ee ee es ee es ee eee ee ee ee ee ee eee ee ee ee eee ee eee ee eee ee 


Oar ie oD ore Ce. SU 


_—— tp GE ee ee ee eee eee eee ee ee ee ee ee ee ee ee eee eee ee ee ee ee ee ee ee ee eee 


Us 3.2 J.974 Ue 2 Oe, 


ae a) Owes rs Oe Ue ou 


ILOhs 





POM Sa 20a tk CASE 4 


TARDE EUR@1 HE. 57H Pew IRXK 


=— ep ee ee eee eee eee eee ee ee eee eee ee ee ee ee es ee ee ee ee ee ee ee ee eee eee ee ee ee ee ee ee ee ee ee ee ee 


CURA EMT NEXT MIYT MUM Cale ay a 
kKS(K) Rone ar REWORK 


—— op 4 eee oe ee ee ee eee ee ee eee ee eee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee es es ee ee ee ee ee ee eee oe 
eR. oe ee ee eee ee ee ee eee ee ee ee eee ee eee ee ae ee ee ee eee ee ee ee ee ee ee eee ee ee ee ee 
=p ae awe epee cee GP a cee 4 eS eee ee Se cee ce ee cep ee ee ee ee ee eee ee eee ee eee ee ee ee eee ee ee oe 
—_ SEP eed a Pee ae ee ees eS SP eee SE Oe 6 ae ee ee ee ee eee eee en ee eee ee eee eee eae ee ee 


—_ ap eee eee, eee ee eee eee SE eee eee eee Oe eee eee See ee eee ee eee ess OS Oe ee eee ae ee SS ee 


Se as. De Os eee Hee 2 


354 Oyen Gens 


—— EE SE eee Ge aee 4 OE ee 6 4 ee A ee A ee eee eee ee ee eee eee eee ee eee ee eee eee eee ee eee 


be I+ a7 o> ier CeU 


Oe Te Me. Dez Use 


Je | i Ue Fd) Jac 0-39 


a a ee ee ee ee ee eee eee ee ee eee ees ee es ee ee ee ee ee eee eee eee eee ee eee ee eee ee ee ee 


G37 Jere Diem Ue FJ 


== <= op ee oe eee oe ee ee ee ee ee ee eee ee ee eee eee ee ee ee ee eee ee ee ee eee ee ee ee ee ee ee ee 


2 cee ce ee ee ee ee ee ee es ee es i ee ee es ee es ee ee ee eee eee ee eee eee ee eee ee ee ee ee 


<< ep com eee es ee coe ee es eee ee ee eee ee ee eee eee eee eee ee eee eee eee ee ee eee oe eee eee oe 


Je S4 Oe) ie Qe2 OR) 


=> ce ee See ee eee ee ee ee es eee ee ee es ee es ee ee ee ee ee ee eee eee ee ee eee ee Oe oe 


ce S33 eo wee ek. Wie oo 


Qe 3? Oe 374 0 2 erecse 


= ee aD eee ce eee OE ee eee ee eee eee eee eee ee eee ee oe ee ee ee ee Oe Oe ee eee eee 


—_— [a ee ee ee es ese ee ee es ee ee ee ee ae ee ee eee ee =_—eame SEP ep Oe ee oe ee eee ee 


106 





COMPUTER JUTOUT  caSE 4 


Peo ferme THe wee we aK 


MUR ES F alae a ae a SPTIMAL 
RFSCK ) RSVR tL Seles) RE WIK K 


a ee ee ee ee ee eee ee ee ee ee eee ee oe eee CR eee ae ee Da ee ee ee ee ee Soe 


enw C ie 7 ee Oe Veo 


cD oe OO ee EP ee ee ee eee ee eee ee ee ee ee ee ee ee ee ee eee eee ee 


Cia 33 0.734 Jez J» 20 


anmenn ae ewe 8 wea aw 2 2 a ewe ae ee Ge ee ee ee ee ee ee ee ae 2eoEe eo PRP aw es Saw =e ae = = 


— —P Pee Pe ewe ees eePe ee ee ee SE eee ee ee Se eee Oe ewer 2S BP te ee em am om 2 2 22 ambamtam=.2. = m= 


a eats eae Je? Jecd 


— PP op eee oe OP ee oe ee ee eee Se eee a ee ee eee ee Pe ee eee Oe 


Ce So er gaet oe eno 


ee ee ee oe EE ee eee ee ee ee eee eee eee eee Eee ae eee eer eee ee SP ee ee aE eee eee ee 


JeGd Yel 43 Je 3 ue DJ 


an Ge a 22 42 aw EE eee 2 SE ae SS aw ean ae 42a 4S SSB ee aa 422) Se Se eee eee Soe oS See ee ies SUC 


Ce. 94 2 i eae O.s reaig 


nm eamapit ap 2242 2 Se awawean 422 enw 22 ae se ee Se Pe eee ee Ss See ees eer ee oe ee eee ee oe eee 


<-> Oo 55 eo Osco 


aw @p a= 224 2 ae 2 ep Sew oe oe oom « aD ee a ee ee ee eee ee ee eee ee ee ee oe 


GeSl Ue e234 Oe 2 Gao 


ew 4a ea 22 ee eae ae ee 2 ae eae ee es ee OE OE 2 ae aewewen SF at 2 Ea CE Ee 2a aw wae Se as awe = —s 


oF 34 Ge ret Wiese U65U 


a oe ee ee ee ee ee ee ee ee ee ee eo ee ee eee ee ee ee ee eee ee ee 


Oe re Geo 20 Je 3 Ue 590 


= =p 4 ee eee ee ce SE ce ee ee ee ee ee ee ee ee ee eee eee eee es ee ee ee ee ae = a 


Cie s OS re Bl age Oe 50 


—_— Pam ete 4 422 ae ae amaet 4 SE oS See ee See SPO eee eee Se ee eee eee eee ee ee ee Cm Se eee el 


Je 3d ae Oe Jes es) 


-_— 2 ep cee Sep SO eee SD ae oe ae oe eee es SP ee ee ee ee eee Oe ee eee oe ee Cm eset eee ee ee ee eee eee ee ee oo 


Ose > Jer w> Ke 2 Oi, 


Je 34 wes CS ALS OG) 


a a a es ee ee es eee ee eo ee ce ee ee ee eo ee ee oe ee ee eee ee 


Jeo 3 Vee nee Oe 3 Oe50 


a oe ee ee See ee eee ee eee ee ee eee ee ee ee ee eee ee ee eee ee eee eee eo 


erode a3 3.7 Ged Wie 3 


Oe8i Oe 75 Qe 3 QVedd 


Ge 3U Oe SS. O.5 Je 50 


107 





Coe oT Es One Tt CAsc 4 


Coote eo He Od Rod JR 


CURRENT NEXT MUN] “ty 4 IPT I Mal 
Hola ont Lo CST RE ACR 


ee ae ee ee Ee ee eee eee ee SP eee ee eee ee eee ee ee ee ee ee ee ee ee ee ee eee ee 


26) Ge Geo OAD 


=p eae SS Pare ean eee ema ean ee ewe Ss TSE wae? Sew eh aes ae eee eS ee eee eee ee ee ee ee ee ee 4 a 


aro 9 Or o53 Veo 0.6) 


Ea eae ee See eae See eee See eee eee eee eer eee ee ele ee eee ee eee ee eee ee eee oe Se See 


1 ¢ a Aare yo . 
Ue ~~ C De 74 oe reno 


Ce >+ OF Ree, ore Ueld 
Geo De =D Oe 3 Q.69 
OG eae eos Ow d Gis 


Cao ae Qe3 Je 30 


—weabamw 2P«e =e ew eT ESE eee ese eee eee aD ee eee eee ee ee ee ee ee See eel oe eee oe eee ee ee 


Diet Sie Ge 2.4 Jed Ores 3 


ap eT SEE EE 2 SeRP ae EEE 2 aE AaB Seen eS en ae TE aS Uae EE ae See eee SSO 


Cec 3 Were ora Je 2) 


Oe a4 Ue ded Ono Oro) 


— <= ae au ee ae oe Ge ae ae eee SE ee Ee ee ee ee ee eee ee ee Oe ae eae ee eee ee ee eee ee es 


0.235 eo 22 Ce 5 Ge SY 





Wmaanrmcarrs SS Ee om eB ew cy 8 OSE Se 2 ee ee ee eee OE ee ee eee = oe 


0.35 a 20S 9.3 SoU 


a ee a eee eee ee ee ee eae Ee Se ae ee 4 eee ee eee ae ae eee eee ee eee ee ee ee ee Ce Se oe 


Ce 34 Je 20 > Ged Je 50 


om ee ee Se ee ee ee ee ee ee eee ee eee ie ee ee eee eee ee eee se ee ee ee eee 


Jeo Ce Dek Digs Cao) 


a ee ee ee eee ce eee ee ee ee eee eee Oe oe eee eee eee eee eee ee oe 


Ce 32 0.397 Je5 UedJ 


= a OP wee ae eee ee eee ee ee eee ee OO eS Oe oe eee eee ee oe eee oe ee ee eee ee SS ee eee ee Se ee Ce 


Oats Ue os Jes i oieal 


mb Ge ee oe ae ee Oe eee ees ep ee eee ee ee ee ee ee eee ee ee ee ee eee ee ee see ee oe 


Ceo & Je 350 Oe. vie oe 


== seep oe Se ee ee ee oe ee ee ee ee eee ee ee ee ee Oe ae eee ee eee ee 


108 





COMPUTER PROGRAM 


ms) DO OOOQOOOOCOCOCOC COCO OCOOCOOCOOOCO COCO OOOO OOOO OCOOCOOOO 
MAN COSI OF D HO AQUA FINVO EK DHOOM OPO THOMNM TINO DMO e1CIM fi Op 00 
QQ0900000 Off tet Sete SI IICQIONICIOQIN AI QUIN gg NM MM CMM 9 OO A ES 
olelel|@lelel eel el Pleo] eee elel e121 e1el ele Fle] e121 ee @1 e121 elel oe] elol lel 21 21ele\elelelelol® 
QO OV VOCS OSSO OOOO SC OCOODDO090009 OOQDADOCOOOOCCOCOOVO0OO0OO 
8 IAA FAP La BE dar AF AT Le RT AF AP LT AM AP A AP LT aT a A AF TL ar a a a ae he” FT la 
De ee ee ee De re Pm ee ee me ee ee ee ee te De ee ee ee re ee 
eee eee ee a SSeS SSS ee Seeaee@ e6eosSsecees¢ 


e 

° O => | 

° = = Om 

® tang regis oo, 

e = Ze OZ 

e ae roy TELLS cl 

e <I =2-==> Li 

° ae ml 4 = 

® Ox = ee yO 4. cs eu 
° a) LL} eee Ls —— Wa OFS wet oie 
® si p= OK OY © Opal fp PT <a Od 
t mt Oo RRO OL DY tau 
° oe) Odgtawn Z2ettot oo 4 > <O 
: = OQ aw mee Fr tft Or we e 
e = NNT D> Ieper UY UW 2 2 wou 
e << YW J CxYtWW < <f —~ US Og =— 
e ad <{ UWIus <f eRe avzZ2ee sP_ wt Wee il 
® > -Y> MAD a mT YW YC Re FY 
° Q @ usw HZOORrW WMO ZerwHtn 
° a reHw COw~aqaonr aul D> o Cao _- 
® © ReuIN GS Ow ce Ow Ww UO>ecaz 
® zz € LL CY wer LU Oe i LS 
e oe ma Oe O2zs = adit WU D> cg WW 
af WY ae 20.2 Wisco co = ad tI GeOwa 
® = Trey C) aao<ti—-C Us us eta 
° ~ Le) Oda CQ te —- cw 
e = << ~tiza Ieqqtrr tet Tw tenw rmnwocd 
° oO VY) SIO eo mye Ie EN ke OC OC TL 
® é < cm CVG Ss 2 Ce> > ss YUL LL UI pe 
° fm = pL <f J I<cf z- Be Lis BD hem 
° z= i (o>) cD >> a me NRL 
e 2 Q ee ULI SLL, LL mnt oo- eI 2Zaqtyj a Oar™ 
e © W = <I. Y) WUE WD Wrew Wh wx 
° 0 < <I ps <I <f Ree SY 2 C£MN> SO+DwW 
* 2 = © Fd pO aIqZZo> oO Yu WwW Ut 
¢ > << <I eZ iL ReROOeU “Tr Daw WWW> B 
© aA = w us NOON WO Oct Wu 
e —e 6 Oe LL AMPS am Zoie ZQ FAWWJS Ep 
e ee a Owe wo WII sees OE 
e a Zz @ Cyn Oe Tres rsYe SF WEIN ANDO 

° oOo a az tf eit Lh Fel pe res OUSLY 
e —~ & us O ZnO Sct OQ— ArRUCOm FENZ 
° Zz ae Oaq-0f- ” tet ie GTOZZ DW el 
: oO Ww ‘5 O aD er OQO0C0OO2Z? CormV%at aerr_ 
e mw | © ot JO us zt S272 =~ D DOZWMesM 
° -— & x Oo aw E Ee CQOO0OFD FHeWOUG Gnato=r 
e ~ WwW OF wf ww IO? OO UW weZZzaw OI OF pant —) ont pm 
° Owe oa ZzorWwo > EE be aa Zdiibmih LOCAdc 
e | mt <t wn | =—O +r I ocoe>r -OrOrFO oO 2a 
e L$ peg b— US pe LU ae OMOMmatoH Aan - © ers > 
e Y Y > AakO TZ «<f Cet CY CVer Cam 

\Y tesa tks, oe) OWZ et Qn OL I re a Le 
® = xr << Us WWI SaoeysrD ONoOoOCcweE 
e Pt = LL $ oes CO oe LS u 4 SO SQ Coed YMWZWS J2eCrvrocd 
e 8) << yy WN TuwwWorw QO 36600 2_~Oa2Oe D> i2a> 7 22— 
? — uw Co — mm bh UL ae <I J- J~ <i >toctZe2z>Zzuormn 
e -—- wd § YW wm W CY Le k= Ss” = 

° qt = OO uu us fw Om <r © 

) ons ££ Tr Taz ZILG= 6 _ a Zz =< 2»x<~-~z 
° EE a ee ee = =e je 2xXZxO - ras 

® LL. Wj tte <= CO. 2 ee <L ey T p 2: 

. rant Y Mm WYO — Sass" a) 

° — 0 On> ww wv <~«x2I350 = 

e > OQ mOM<_™~ O oe) 

e LJ a cl Fa) 

@ co =) ut 

e aa fae C 

e 

@ 

@ 

WO 


OUI OY QO OOD OD DO QU U0 DV OO UOOD OUUOY OO VOU OOO YD DOU DO VOVOUUUHO 


Oe 





“ 


EVEL OF THE STATE 


Mo OODDOVOO9O ODO 9000 OVO VO 9OO COC OO DOC OOCOCoOO9OoOoO90000000So 
DOMINO SFO0 MD HOM QyCAT WOM DDO emN Pin OM DHS sayM.¢ NOP DOMOMOIM TINO 
AP LAUD LOLI LOLI UV UN £50 5 0 080 6980 0 SO Be Pe fee Pe Pe PHP DW ODO MDDDNDAOGARHANHM 
DOOQVDOVND SOO OO QDDVOAWOTOO OSC OOVTODWOAOGOODGVOOO OOOO COOOCoO00 0000 
OODWDVAO OV OVO OOO OVO OOOO OOOO SOCOTVOOOOOVVOVOOO COC0O00000 
ee eS an e282 222272222 22222 ZIZZO Z 
De ee ee Pe ee ee ee ee 
eC) OO CE OD OOD8 QD QNOOO OD OOOO GOOG AOA OOOOGOOOOQAOAaaAaoG 


e AT THE 
0 T 
V 

L 


EVEL OF THE CONTROL 
E 
H 


CRETE TIME 
p 
E 
A 


—! Om %/) 


OnesnSarw «1 
SU FU Se od 
rg et ey eg §=— CL) 
ed el lO 
tel <T t+ S| 
Or oOruUd 2s 
Pond emer be Se Poy fre OO 
Oo fF and UOae 
ai 0 ah EE be 8 

tet ed po ed 2! 
oo azn 
Dee <TD Te pay OSE LY 
Lhd ee ht] ot fy fam ene teeny) 
O¢OCO—-eEC 
zwatoastaszzvoaast 
best > eet > seg a 


YOUU O00 DUO OU VUDUOU UO UD VOQV YOUU OQYVQYVQOOUVUOYUUOV UO 


CATA CARDS 


LI LL oo 
COle— eLotO 
aS ape Od ek | 
N<fO— j 
=e OQwaw 
rOZWOO « 
Pa oh ee oD) 
Iw we i 
> aed OL et = LL Ot 
On mm 
qUn<t uy 
LL) vat weed OS cl = 
em yO >= 
—Iqt ww D 
re >See oe 
WY) LO7O 
— Jk (90) 2. 
uw OO me 
Ort CwWX<s 
PLL wee Tm <L ey 
C7Z7Q F= oe 
WoO 
MOZNAIO=— 
= }OzZ= 
™UOU <I e 
a Oe = uy 
Oe 
eg MWe 
Ly ot J ef 
tat >-— © 
VZO ew ZiL 
—w el ADO 
VLE=ZZen 


TON LABEL IN COLUMNS 1-72 


Cw) I Om 
WUOMPZz=wen-e 


SECQND: 


AND FUNCTIONS REQUIRED 


SUBROUTINES 


sEILO, 


(5 
rll 
LLS pom bame oe 
po) ae) 
Oawn LL 
aes <{ 
Orn = =TrYN 
Zw eRe 
<—T LL VY) 
a>. WNrO 
ee CO WMYY 
Ot ¥~oO 
emt eT COLL 


US OO Auld = Ly 
xertct HOM 
—E CN 
2 CWoar—y 
QwiWoOAa<|d <f 
WIMNNDNWe 
re Hee Cw 
—- Cc Ke> 
ZrIWOmwW W 
Secrttovnn 
Wot RL ud 
zo Zu 
> a Pade 1 ns 
QUIMeQoa ze 


eee @voeeeevaosvnsesesesveeeeveeeeeeeveveeeveeeeveeeee8eveeeeeeeeeeeeeeevee eevee eev eevee ee 2 eee @ 


L144 y6Xy9ABe///) 


) 
’ 


OOM 





fm OS VOIOOCO CEC Oooo Ooo awmnoOoCCeoooCooooCeCcoCoooooooOo oOo ooe 
DOM ODS stu FMOOM CS POSMAIM SON OF DNS HQ M™POOH DKNOmINQATN OF OCeOnN 
DAN OOOO O00 00 Ce et st 8 1 II IOI N CUIN AIOIQIONIQIOICN MOM EHF) COOMA 
OD OS eed eid neg 04 ced ed gee ge) ed rd red ed med crm) eed) a ed gern) ed ee reed) 0) ed reed greed ge) ed ge Fd cred ed vd ge) Pd ped moe 
PIOOOCOOOQOCCOCOCeO Coe ooceecoooCooosaooseceoooCoeoGooocoOooeooO 
Be a Pee OPPO AO COCO EOC CeO eee EAeeeee 2e2 CeZeZaeeZe 
Spe ee ee ae ee ee ee ee De ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee 
DADBAQQANDOD QOADOD O08 QOO COOH QAD A QOOGCOQOKOAOQOQAOOOQAQ00 


Fouy) 


LL) -” = 

OQ = ~ ew 
© <I a at a) 
= = od —~A 
= © -— 
om i) o = | *< 
<x<- OC mW - = 
ax ¢ Ne ~ 
*<ft = ow ny 
2D OC 2x ~o 
- & = & ai b= 
ox J me my) 
<et Teg rc) 
Oe SNe WY “NO 
rw a ae Oe 
>a fe we - & 
<m OS wo Uj >< 


De YEO WEO 
OSS RKOOwtk™ Cee lL 
ax“ OMiLs 4 & eos ® 

ata “NO PO > ow em dl Oe 

e OU se Mew wo FZ 

<I <ete US © ©] bore © 
marie aeDNX<Ale =e © oi 

— IML <T OM eww XK Ct 
M —O aelinme “qomnnonwo 
e ete TNs & ML TOAO ML 
Oam @ OX © @& RN © BO Ow ~ 
SHILA © xXO>xe ed | OX 
Le LO Ne Ope Neto TKO 
NeapstoaMeqme JNA mat Mm 
et ae nes SP gt SO ee TL tee eee SP ee ed Oe 
hep bape Om RO LS pe ee Le 
ele <f x c(a( per e(dicte <f 
a be a oe i ee oe 
ae Tama Male de ee Gas an A ta ae A 
OD0O0D0 pm0Ic0sOoOocoeosro 
Kee bag ee hn NL ay OMe tL LL LL OL. 
ot ets and 


COon oe OO mW 
TaN AMCs 


DATA INPUT 


© 
TUAW OULU 


— mom KOO 
~~ a OS om MCUCUCN 
(BS apemiQnies © & © & 
mt oe & ww MOO00OO 
eoklattallallat Ahem)! 
wee were er Lay i 9 
AOoOookr-e-<t 
frm XT cl SL st Peet peeey rot 
<T LY Utd LL or Of Oo Oo LL 
OCetxrxrseessaes 


aed: 


UC XMENgXMAXe Ne NXe JS 9X eQUANTX ) 


LO et ad on) ae er SO) 
Wit+Ux<a~— 
ony >t w OU 
O92 YZ RR 0 

— 1 JY) nd Se wee 
OWOrm<t 1 Om ™ uw 
OO ae I) Se 
oo <a 
tt ttacg@ Jie 
YWDFzO Os 


51) 


ITIALIZE COST OF THE LAST STAGE 
MAIN LCOP 


IN 


© 
co YOUU 





Pes SS) So OO OOO OOOO OOO OOOO OOO OO OC OO OOOO SDO BOO OOOOoOoOs0ooO 
PVP DO PP DDS eg QM FLD OPED MH Meet MO PUNO DWH OM QUO SIN OH DODO FACE ONO 
ade Ag ge SR SALA LALA LA OA LUNN ODO OW DOOR EEE EEE EP DpDDDDMNMDDOMOM 
etd med md ged ge ges ad md eed Oe gee nd gm ed mer] omen) ee nd ound omen) gud gre] rane) goed) eed tt me] em om) reece unl) geet geed] geen PA me gaan) OA cpr ome ge 
Pe i OOO O00 90 COO O COO OO SOO OOOO OCOOCOOOCOOOOoOCVOOOoeo00o00O 
a i ee OO PO EO COL CMP EEC EUEECBEBEULEMECLLLECECOCBECCECOCEE 
De ee ee ee em ee ew em ee ee ee ee ee ee ee ee ee ee ee ee ee > 
DOOD DAHA DODDOO DD AOOCODADOOQO OO AD OOOO ODACOOSCOOOQORQQAD00 


@ | 
e e ae 
e @ a 
@ e LL 
e B Cc 
e e z=. 
Pol @ e —¢ 
el e e & 
i e e pan] 
ad e e + 
~ ® e ww 
ae @ ® eo 
- e e oe 
as ® e fj) 
on) e ® = 
cad = e e ) 
— ee @ e in © 
o~ am e @ i) 
Ss tmesnd e e CO + 
- 4 o e Re oO 
eis = a 8 e <I 
— <— < om @ Cpe ot 
a. Zz > ~ (9 * ~™ C5 LL 
> 4 < ea) Co i eos rag Hf 
p= = ~ & om Le eae (_) me 
= Co yd (_) ao Cte @ eS) palate 
<{ = QO a YY) om @ C2 =) 
= = >< md © er & — tty <P OW 
| _— To i) e> © ww = of SIO) 
o =5 << ~ eae © eee ee tiie 
a <i =) = ene of <{ <InuOm 
& “J le eo. enn] © a > — ® ~ 
= “i Oem © eo & SO OL el 
> anton} OO os ex 6 ws ey “uJ + 
S =—_* 2PO oe eee Ly ~ op 
ae ~ 7 eZ eed >< eH e* oc <~=—O—~ 
oO Oe x Meo ~ eee © 8 — Oz 
- NZ Ir <~O x< eze D = 2e O 
om ~ & = §*OD eUW-e ~~ emmy ® & W) tyne 
i O>< Det <e —) ose = DULL 
~ ws 2ax<O “79 = ene © i— % Det 
~ he >< — MHZ Se O exe mt < — 00] 
—_ — o =D ee Zw) — eo @ = — gard OD ome 
_ (= D Steet ed et fl ae) ee#* O o. — 
ie 8) TZ bout: we fF oS TD — iT oe) —J oll + 
¥ a ae OY OVU ws OT — ene QA We Woy 
aad OM DeWwW— WOR, = on om §§ CE COO Oc ww - 
— om OR = ene ee JO mz *2O W 
aed a> eteftOntw *wOorR< nme ev*® wom twat +s JO OO 
mt UE ONMD RE TORE FIO OS em © © <meta 
See ot OM OOOO Zee 1 8 em © OF SW FARRER 
m~elOotas w =F eO™wOOUWG ower NN 8 peta 
medal a yt STOTT w CO Or ex q_JjOwWos+ —_ 
CAC Ye JOO Ld we ad we MOND LN ed? CO ™® A2DWUNOY + 
eo Si mw Ocdtef Lei LbLOnOaLastuiu & e e te UrCOOuU a = 
OO Sei MQM men OD Hh OO ew? ££ OOlzn—~!1-Cfa ~ 
eer OS I~ WY) ZoOoae om © W wl XM REY WWW Oo 
P WW Set LOO Os a LW ese QO yo mcw FOO 
Mh age) IO Ly tb ee © NK WWI OZM ah 
pany md we ey ema? 2 KNOY UKHO NOT 
lero WL O ct a cme zane COO+a 
~SFzO m4¢_) Fp GO. er © = et NW tek ee OD 
zoo ome 2S = SWOT ™ If Or 
Ot = e™> eo OF ft Wwe ehL sO title 
Owmuw one VY CF DAO™L=| OV Q0O%W 
td] 
i) 
2 
2 
® 


© 
MH) OW OF YO WO L\ emt 


2 





e, 


CEO COO) 8 COSC O'S 010 OO O'S OO: 0:6 OOS O00 C1010: S C1010 20 10)'0 0 0:010, 0 0 0:0 0:0 000 ee 0 0 ee Oe et eee) 


fINYNO SS SWISN INS XVWSS NEWS NANVNO JNIINOVENS 


@eeoe0eeeoe3454u@?;8eseeeeseee1eeseeseeoeeeeeeeseeseeeoeoeeeteeseeeoeeseeaeeenreeoeseenwreoseeensneoeses ese eeeee 6 Ah a a ic | 


J 


u 
GS T104N xO 


e€*x*N=9 >A TdWVXS 
° 33H INdNI 3¢G OL SVH 9 NOTLINNA SHL 


WU YO 


@eeeseeeoeoeeee820e2802028080 0 Oe Oe O48 OCS Oe 8). 0 OOO OO Oe Oi ee Fee OOOO OOO eee a eae ok Mee eee ae) 


(N)D NOILINAS 


eee eevenevseecneeveevoee ee Cee eneveveeeeevneeeveeoea eevee ee eoeevoeaeese ee eeeaeeenenen ened 8280800800 mk 6 


a, 


0° 000 Gx ( 


e**X=NH 257 


WOUWOVUYD 


Coeoerceresreee reece eereoreeeeeeeeerereEDEeSEeeeeeeeESD EER ee ESET EH EES EH ESE SLES BO 0%) 


(X) NH NOFLONAS 


Cee 270 0 0.0 O:0 016 26 616.68 O16 0:4 0 CO 018 e 0 AS:6: C6 0 0:08 CO 808.0 OO 0 0 Oe Oe roe OO OO ee ee) 


2 


NxLWo~nd + X*x( 1739 


GxV 
38 OL SJAVH S 


Oa 
sd 
= 7 


+ 
N 


i*e@*xyvwn'n INN XVWX*NIWX/YULSNOD/NOWNOD 
CSCO EHOCHT OSHC EET THESE EH SCH SH SCOEHHHHEOHEE SES HOS SFOS 02 010 00 000 20/05 88 +0 2 > 2S 


COVTSATIONGE WS TdxXPN SX ALVLS ANTLNOVENS 


SC SSCHHSSHSHESSEHSESHEHSSHESSEHEEHSEH HH HSEHHEHEHEHEHSHSHEH OCH SHSHH OSHS SESH EHSHSHSEHESOHS LORS ESN 


DADDQADAG GOAQSADQDDAAGSHADADARDCAGQAADAQAANAADAAaAdnm aa CAaao 
WO OWOUUU 


mao IOS COCO COCO OOO SOC OC OCOOCOCOOCOCOCOOC COC COCO OCOOCOCeocoO 
ed NNT LN OF 00 SS eg 0 LO Pe OO Qn met QylASF LN OP DA OMA MGM FINO ODO HOM 4+ OF 
DOADXADDUMH CSOD O OOOO O OC O68 lett td MRI AIC QININ AIO Qu Qin Mm Meg MMM 
rr gd ed ed md rd teed md CUI CQUION OY CQUIRQUIAQU IOQUIQIOY OI QE QIAO AI OI NOU CAI QU CAIQIQIN QuMialAiauMiaqiN qaqa 
DODO COCO TOSCO DQOCDOSCGOV DOSS CSCS GOVO SOSOVDOCCOCOCOCSOOOOO 
te a ee 2 2 222 222 20202020 2022 202222222222222227 
De Da De De De De De ee ee ee ee ee me De ee ee em De Dm De ee ee ee ee Pe De ee ee ee Oe 


Lis 





114 





10. 


BIBLIOGRAPHY 


mupecme B. Bishop (1975), Introduction to Discrete Linear 


Control — Theory and Application. Academic Press, New 
MOT . 


James Guy Bohannan, A Reliability Model for Stored 


Items Requiring Rework. M.S. Thesis. U.S. Naval 
Fesegracuate Senool, Monterey, California. 1977. 


Donald E. Kirk (1970), Optimal Control Theory — An 


introduction. Prentice-Hall, Inc. Englewood Cliffs, 
New Jersey. 


bemgamim ©. Kuo {19/77), Digital Control ‘Systems. 
SRL Publishing Company, Champaign, Illinois, 61820. 


Richard E. Bellman and Stuart E. Dreyfus, Applied 
EVeomee sProgramming (1962). Princeton University Press, 
Princeton, New Jersey. 


MOOG sR. V-, and Craig, A-T:,) Introduction to Mathemacical 
Statistics, The MacMillan Company, 1959. 


Feller, W., An Introduction to Probability and Its 
Applications, John Wiley & Sons, 1957. 


Glenn F. Lindsay, Reliability in Stored Ordnance: 
Selection of Items for Training Expenditure, Naval 
Postgraduate School report to OP-964. 


Donald E. Kirk, Stochastic Control, Notes Naval 
Postgraduate School. 


Ogata, K., Modern Control Engineering, Prentice-Hall, 
iicwe sng lewood Cliffs, N.J. 


a Hales: 





Pat Drs TRI BULLION blot 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Bribery, Code 0142 
Naval Postgraduate School 
Monterey, California 93940 


Associate Professor Glenn F. Lindsay, 
Code 55Ls 

Department of Operations Research 

Naval Postgraduate School 

Monterey, California 93940 


LT. Ellen Franz Roland, Code 55Ro 
Department of Operations Research 
Naval Postgraduate School 
Monterey, California 93940 


CDR Paulo Antonio Ferreira 

aos cuidados de Antonio J. Ferreira 
Ree brs OGtarilo Kelly 20. apto. 403 
Tijuca Rio de Janeiro 20000 

BRASIL, SOUTH AMERICA 


iMiele 


No. 


Copies 


2 














¢ & 


thesF2647 
Reliability contro! model for stored 
Ty) en Te id 
| | 


| 


yi 


AE iit 
i | iH Hil 
TOUTE 
3 2768 001 01500 1 
DUDLEY KNOX LIBRARY 


ite 


| 





